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MILITARY PROBLEM 

The combination of shrinking financial resources and the prospects of a smaller, 
all-volunteer Army will increase both the demands made on Army personnel and the 
importance of the individual soldier. There will be a greater need for more effective and 
efficient training, adequate to the task of providing an increasing r.umber of complex 
skills to wide'y differing students, while using fewer skilled instructors. 

The most promising approach available to meet these new training demands is 
computer-administered instruction (CAI), if it is developed as a comprehensive, total 
system. 

The goal of Project IMPACT-CATALIST is to provide the Army with an effective, 
efficient, and economical CAI system in a total system framework. To be effective, the 
system should maximize the achievement of the students and the instructors to a greater 
extent than is possible in the traditional classroom; to be efficient, it should provide 
maximum productivity per unit time on the part of instructors, administrators, and 
students; and to be economical, the cost and resources must not exceed those of a 
comparable effective non-CAI instructional system. 



DEVELOPMENT PROBLEM AND APPROACH 

The objective of the CHARGE Interactive Graphics System effort was to expand the 
visual link between the CAI system and the student. It was believed that the overall 
effectiveness of the training process could be increased if the CAI System were provided 
with the ability to logically construct and present to the student graphic images of the 
most general nature. 

Accordingly, at the outset the CHARGE development effort considered graphic 
display of dynamic 3-dimensional real and abstract worlds, A system able to handle 
such problems can easily handle subset problems such as display of alphanumeric text, 
line drawings, and 2-dimensional graphics. 




RESULTS 

A powerful interactive graphics system, called CHARGE, has been designed. The 
CHARGE system has three components: 

(1) The host central computing system that manages the CAI system and 
maintains a symbolic representation of the 3-dimensional world under consideration. 

(2) An image generator functional unit, time-shareable among the CAI users, 
that translates the 3-dimensional world description and student responses into an encoded, 
colored, pictorial image. 

(3) A graphics terminal that buffers the encoded image at low cost and 
displays it to the student as a colored pictorial image. 

A detailed description of the 3-dimensional world description and image generator 
function is provided in this report. A description of the CHARGE terminal is given in 
HumRRO Technical Report 74-26, CHARGE Interactive Graphics System Terminal 
Theory of Operation, December 1974. 
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PREFACE 



This report describes the results of the development by the Human Resources 
Research Organization of the CHARGE Interactive Graphics System image generator 
and author language support. The image generator is described from a theoretical and 
detailed functional processing point of view. 

The work was begun at HumRRO Division No. 1 (Systems Operations), Alexandria, 
Virginia, under Project IMPACT in 1971 and has continued under Work Unit CATALIST. 
Dr. J. Daniel Lyons is Director of the Division (now the HumRRO Eastern Division), 
arid Dr. Robert J. Seidel is Program Director, Instructional Technology Group. 

Members of the CHARGE System R&D Team were Dr. Ronald J. Swallow, 
Principal Investigator, Mr. Roger L. Gunwaldsen, and Mr. William G. Underbill. 

Detailed description of the CHARGE terminal is given in the report: CHARGE 
Interactive Graphics System Terminal: Theory of Operation, Technical Report 74-26, 
December 1974. 

The work described in this report was conducted for the Department of the Army 
under Contract 19-73-C-0004. Computer-assisted instruction research is conducted under 
Army Project 2Q763731A734. 

Meredith P. Crawford 
President 

Human Resources Research Organization 
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INTRODUCTION TO THE CHARGE SYSTEM 



BACKGROUND 

The objective of Project IMPACT-CATALIST is to evolve a series of prototype 
systems of computer-administered instruction (CAI) in order to produce a total CAI 
system that is effective, efficient, and cost-effective for operational use in a training/ 
instructional environment. The total prototype sy tem includes four main components: 

(1) Hardware-Hhe computer, student stations, and related equipment. 

(2) Software— the computer programing systems that control operation of 
the hardware. 

(3) Courses of Instruction-the actual content and logic of courses administered 
by the computer. 

(4) Instructional Decision Models (IDMs)-the rules and strategies by which 
specific course content is provided to an individual student. 

As part of the Work Units IMPACT and CATALIST, HumRRO conducted a system 
design study which led to the development of an interactive graphics system called 
CHARGE (Color Halftone Area Graphics Environment). Designed initially for use in 
instructional systems, the CHARGE system can provide considerable and versatile graphic 
displays for many applications (including equipment design, troubleshooting, etc.). While 
being cost-competitive, the CHARGE system can achieve color, gray level, three- 
dimensional perspectives, and rapid updating. Perspective views can be simply, rapidly, 
and economically generated from a (3-dimensional) xyz-coordinate, dynamically changing, 
"real world" description. 

SYSTEM OVERVIEW 

The CHARGE system has three major components: a hcst central computing 
system, an image generator, and a set of user display terminals. A detailed description of 
the user display terminals is given in the IMPACT-CATALIST report: CHARGE Inter- 
active Graphics System Terminal Theory of Operation^ These three components are 
configured as in Figure 1. 

'Ronald J. Swallow, CHARGE Inieractwe Graphics System Terminal Theory of Operation, 
HumRRO Technical Report 74-26, December 1974. 
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Figure 7 - The CHARGE System 



The host central computing system interacts with each user as he creates his 
3-dimensional "world," a description of a 3-dimensiona region of space. The host does 
this by maintaining and manipulating two separate files that serve as the data base for the 
CHARGE system: 

• The first file is the "atom" file, a low-level data base that defines a library of 
standard 3-dimensional shapes or building blocks such as a unit cube, unit 
sphere, and so forth. These atoms are generally not modified by the user; in 
fact, the user does not even have to hnow what data are used to represent an 
atom, he merely has to know each atom's standard orientation in space. 

• The second file is the "object" file, a tree-structured file by which the user 
can create and name complex objects by defining them in terms of pre- 
viously defined objects and atoms. The user does this by use of symbolic 
commands such as: "Create object A by adding objects B and C. Then create 
object D by adding objects A, B, and A, this time moved up 3 meters in y in 
the xyz space." The linkages between levels of object definition are "trans- 
formations" (such as, "move 3 meters in y"). The lowest level definition in 
any object is always a transformed atom. The object file structure for object 
D in the above example is given in Figure 2. 
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Figure 2- Object File Structure 



Atoms in B 



It should be noted that the data structures in the host central computing system do 
not represent the graphic view of the world; rather, they represent a data base from 
which an infinite number of graphical perspectives can be made merely by the command 
"Put me at (x,y^) in this world and show me what I see.'* In fact, the host can put 
several users into the same world and generate views for each using the same data base. 
In this aspect, the host system is far more efficient in its file storage than it would be if 
it had to store **canned** (i.e., pre-computed) perspectives. Furthermore, the CHARGE 
system has the capability of generating a dynamical representation of this 3-dimensional 
world by giving the transformations jnd/or observer time varying parameters. The cost 
and performance of the host is shared by the N users. 

The process of generating the perspective begins at the time the user gives the 
**show** commeuid to the host. The host then initiates this process by first **compiling" 
the user*s world. The compilation process accomplishes two tasks: 

(1) Each atom that is symbolically referred to in the object file has its atom 
data retrieved from the atom file. 

(2) Each transformation string linking an atom to the world is replaced with an 
equivalent composite transformation. 
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These two sets of data now define the user's symbolic world. There is one 
composite transformation for each occurrence of an atom in the user's symbolic descrip- 
tion. The amount of actual atom data is related to the number of different atoms 
symbolically referred to by the user. For example, if the user's symbolic object file had 
50 atoms, >,ach one a Cube, the compiler would produce 50 composite transformations 
but the only atom data that would be retrieved is that of the atom Cube. These data 
along with the user's viewing orientation in his G-dimcnsional world are then sent to the 
CHARGE image generator for perspective computation, thereby freeing the host for 
another user. 

The top half of Figure 3 dcpiv^u. each user's i\o\\ through the host computing system 
as he creates and views a 3-dimcnsiunal world. A typical user will spend most of his time 
in a **think" mode causing the ho^t to loop through all of the **No" branches each time 
the user is polled. The tup three boxes arc entered infrequently compared to think time. 
When they are entered, they art tvpicdl interactive user requests requiring very small 
amounts of the host's resources. 

The greatest load to the host computer is that of compiling a user's world, but this 
too is a low-frequency event. Furthtrmore, the demands on the host are still quite small. 
A world with 100 atoms, each wii an average linkage depth of 4 levels would require 
approximately 100x4x27 - l.lxlC* multiplications to produce the 100 composite trans- 
formations. If the host has a multiplication tiuie of Sjuscc/multiplication, this would be 
33 msec, of computation, a fraction of the real time required to retrieve the user's atom 
and object files from disk storage. The computational load of computing the user's image 
is off-loaded onto the CHARGE image generator. 

The CHARGE image generator is a processor located between the host central 
computing system, and the N ustr tenninals. Its function is to translate each user's 
3-dimensional world/viewer data into a high resolution encoded image that can be 
buffered locally at low cost at the u.scr\. terminal. Like the host, the image generator's 
cost and resources are shared by thi' s>stcm\s N users. At the present time, the image 
generator function is being emulat<Mi in softN\arfs thereby allowing for experimentation in 
new algorithms before a hardware design is finalized. 

The CHARGP> image genrrator is reallv a t\\ onstage processor. The compiled world 
data from the host enters i\v projei tion proci^ssor which projects each atom into a 
"perspective domain** and hinits thi n suKing pcrspi-t tiv<- domain data to those atom 
surfaces that would be visible- to thf user on his monitor if no other atom surfaces were 
present. The exact form of this piTspei tivc domain atom clata is chosen to maximize the 
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rate at which the second stage processor, the visible surface processor, . able to 
determine which surface is actually visible at each {x,y) coordinate on the user'^ .nonitor. 
The visible surface processor also encodes the resultant image into a highly data- 
compressed format that can be efficiently transmitted to and economically buffered at 
the user's CHARGE terminal. 

Once the encoded image has been transmitted to the user, the image generator is 
free to service another user. The two stages of the image generator are independent and if 
a second buffer for y-sorted perspective domain atom data is added, the projection 
processor can fill one y-sorted buffer with one user s data while the visible surface 
processor empties the other y-sorted buffer containing a second user's data. 

Figure 4 d^^monstrates the operations of the image generator on an extremely simple 
case. The user has symbolically created a world consisting of just two atoms: a cube and 
a triangular plane. The host sends to the image generator the following data: 

(1) The user's orientation 

(2) Atom data for the standard cube 

(3) A transformation orienting the standard cube in the 3-dimensional world 

(4) Atom data for the standard triangular plane 

(5) A transformation orienting the standard triangular plane in the 
3-dimensional world. 

The projection processor first projects the two atoms onto the user's "viewing 
window," the rectangular area corresponding to the user's monitor. Back surfaces are 
"clipped away." If any atom surfaces had projected outside of the viewing window 
rectangle, they too would have been clipped away. The resulting perspective domain data 
consisting of surfaces and edges is shown in Figure 5. Note that no "front surface" 
decisions have yet been made between the surface of the cube and the surface of the 
triangular plane. 

The visible surface processor resolves conflicts between these overlapping surfaces 
and encodes the image into the edge form of Figure 6. Note that as conflicts between 
overlapping surfaces are resolved, intersection edges are automatically created. 

Figure 6 also demonstrates the efficient data format at the CHARGE terminal. Each 
"visible edge" in Figure 6 defines a color and a brightness to its right side, the edges 
themselves are ordered in x such that for any .stanline (y position) all of the edges to the 
right of Edge i have edge labels greater than i. Edge 1 is a special "marker" edge that 
flags to the CHARGE terminal decoder unit that there are data present. The image of 
Figure 6 contains 20 edges requiring that 20 76 = 1520 bits be buffered at the CHARGE 
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terminal to produce a sharp, colored video image on a monitor with resolutions 1600 
units horizontally by 1200 scanlines vertically (4:3 aspect ratio). If the same image were 
to be buffered as a video image with 9 bits of brightness and 12 bits of color information 
for each point, 20.7x10'^ bits would need to be buffered. 

The CHARGE image generator and terminal also have the ability to "smooth-shade*' 
curved surfaces; that is, curved surfaces are approximated by planar facets and the 
resulting edges have bits set instructing the CHARGE terminal decoder unit to linearly 
interpolate the brightness between successive edges, thereby giving them the appearance 
of a continuous, smooth, "round" surface with no edges within it. Thus, round surfaces 
such as spheres and cylinders appear round on the terminal. 
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Section 1 



CHARGE IMAGE GENERATOR 



OVERVIEW 

The image generator may be thought of as a very fast, highly intelligent data 
translator located in the output data path from the central computing system to the 
display ternunals. Input data from the central computing system are: 

(1) World description: A complex, symbolic description of a region of 
finite 3-dimensional space. 

(2) Viewing parameters: Definition of an observer with respect to this world. 
Output data to the display terminal comprise a highly data-compressed representation 

of this world as seen by an observer with the specified viewing parameters. If the world 
description and/or the viewing parameters contain time-varying elements, the output 
becomes a real time sequence of these encoded images— that is, a movie of a changing 
world as seen by a moving observer. 

The image generator is a two-step processor and operates as in Figure 7 on a 
symbolic world-observer pair at time instant tj. The projection processor uses the observer's 
viewing parameters to map the symbolic 3-dimensional world description onto the 
observer's 2-dimensional display screen. Since projection maps each item in the world 
independently of the others, the projected world data must be further processed by the 
visible surface processor to determine which of the projected items are visible and to 
encode the resulting video image into a data-compressed form to be sent to the dis- 
play terminal. 
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Figure 7 - Image Generator Processing Steps 
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This report on the image generator will contain three parts: 

(1) Specification of a world description and viewing parameters 

(2) Projection processor 

(3) Visible surface processor 

WORLD DESCRIPTION AND VIEWING PARAMETERS 
Introduction 

This section provides (a) a description of the 3-dimensional space, (b) the method of 
rep/esenting an observer within this space, (c) the method of representing light sources 
within this space, and (d) the method of representing 3-dimensionaI, automatically shaded 
objects within this space. 

S DImensional Space r3 

All positions X = (x,y,x)T = xT+ yT+ Zk in the world are defined with respect to the 
righthanded rectangular coordinate system of Figure 8 where unit directional vectors 
i = (1,0,0)T J= (0,1,0)T and k = (0,0,1)T satisfy the vector products: 



ixi = 


0 


jxi = - k 


k xT 


= j 


Tx] = 


k 


JxJ = 0 


k X j 


= - i 


i X k = 


- J 


j X k = i 


kxk 


= 0 




The point 0 » (0,0,0^ * Ctt 4 Oj -f Ok is the origin of the world. The 3-dimensional 
space is defined to be finite and will be denoted by R^. 

World Description: Observer 

The observer is completely specified by the set of viewing parameters 
(Xyy,Y^,Z^,X'j;,Y'j;,Z'j;,Xy,Yy,Zy,W'x,W'Y,MAG,BQ) which define an eye, 
a viewing window, and a point of observation. These parameters are pictorially repre- 
sented in Figure S.Y =» (Xy,Yy,Zy)'^ defines the point of observation with respect to 
0 » (0,0,0)'^, the origin of R^. This is the point toward which the observer directs his 
vision. The observer sees into R^ only through a viewing window. This viewing window 
has a center located at"W = (Xyfj^Y^^Zyf)^ with respect to the origin of R^, and lies 
in a plane orthogonal to the vector V - W. 

W is the origin of a second rectangular axis coordinate system with axes x',y', and 
z' and unit directional vectors!',]', andlc' defined so that: 

T'xT' = 0 7'xT' = - k' 

TxT = k' T'xJ' =0 

T' xk' = -T' J'xk' =T 



k'xi' = j' 

k'xj' = -T' 

k' X ic' = 0 



' llV'-Wll 

- - T'.T="o 

With respect to this primed coordinate system, the window size is W'^ horizontally 
by W'y vertically, and the y and y' axes lie in a plane. The observer's eye is located at 
E' = (X'j;,Y'j;,Z'j>)'^ with respect to the primed axes. If X'j> = Y'g = 0, the eye looks 
at V = (Xy,Yy,Zy)'^ directly through W = (Xyf.Yyj.Zyj)^, the center of the viewing 
window. The observer looks through the window with a telescope with magnification 
factor MAG > 0 and coefficient of optical transmission G [0,1] . The relationships 
between the observer, the viewing window, and the point of observation are pictured 
in Figure 9. 

Specification of the observer defines his **cone of vision*' into the world. Only 
those objects within this cone are potentially visible. One further specification is made; 
Only those objects that lie on the far side of the viewing window will be allowed to 
be visible. This reduces the observer's vision to the ^'pyramid of vision" (Figure 10). 
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Figures - Observer Viewing Orientation in r3 



World Description: Light Sources 

Point light sources in r3 are defined by the triplet Lj = (^i,Bi>ro-) where Xj repre- 
sents the location of the i*^ point light source, B: > 0 represents its intensity, and r^ 

OJ 

defines its brightness fall-off ratio. The source radiates its light evenly from its 
point location. 

The sun light source is defined by the pair (Sx,Bg) where S represents the light 
direction and Bg > 0 represents the intensity. This sun light source provides parallel 
light rays throughout R^. 

World Description: Object Definition Requirements 

The task of mathematically representing a complex object and all of its physical 
characteristics in may be quite formidable. This problem is simplified if this 
representation is limited to those parameters required to dynamically position the 
object in R^ and to view it from any orientation in R^. 
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Figure 10-- Pyramid of Vision 



Consider an arbitrary object. The object is said to be ^onver only if the region 
of space S included within the object can be compressed as: 

S = {(x,y,z) lfi(x,y,z)< 0,1 = l,...,n) [i.i] 
whei^e fj, . • • ^ fn are real-valued functions such that for i = 1, . . . , n: 

fj(x,y,z) i = O^in = i fl.21 
I < 0,if j i= i 

define the object's n boundary conditions or surfaces. If an object is not convex, set S is 
. the Union of some m convex sets of the form of [1.1] . S is, in general, an infinite set 
regardless of whether or not it is a convex set. 
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If the object is to have a dynamical description in R^, set S must be given a mass 
density function dM such that: 

M= ;dm [1.31 
S 



\^ere M is the total mass of set S. If the object is rigid, a point mass representing the 
total mass M of the set S and I 
The center of mass is given by: 



total mass M of the set S and located at the center of mass G = (Xg,yg,Zg)'^ is sufficient. 



XgM = /xdM 
S 

[1.41 

■ygM = /ydM 
S 

ZgM ' = / zdM 
S 

If the object is homogeneous (i.e., dM = constant), the center of mass G coincides with 
the centroid C = ix^,y^,z^)^ of set S given by: 

V = /dV 
S 



= /xdV 
S 

y.V- ;ydV 

s 

= /zdV 
S 

where dV = dxdydz is the differential volume element of R^. 

If an object is not itself rigid, but is made up of n rigid subobjects, then 



[1.51 



m m 
= 2 Mj V = 2 ^ 

i=l i=l 



m m 
= 2 Xg. M; XpV = 2 x,.. 

i = l ^' ' i = l'' 



m 



ygM = 2 yg. Mj m 

i = l YcV = 2 ycj Vi 



[1.61 



i = l 



m 



ZgM = S Zg. Mi m 

i = l ' 



i = l ZcV = 2 Zc- Vi 
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Since the only parts of an object which are potentially visible to an observer are its 
surfaces, the object at any instant of time may be adequately represented for display 
purposes by an appropriate description of its surfaces-Hhe n boundary constraints 
of [l,l]--at that instant of time rather than the infinite set S itself. Color, opaqueness, 
and optical reflectivity are natural properties of surfaces and will be included within the 
surface defmition. 

It should be further noted that for homogeneous objects these n boundary constraint's 
define the limits of integration of [1.3] and [1.4] and, hence, define the center of mass 
of the object and the mass of the object to within a multiplicative constant. 

In summary, an object definition requires a defined mass, a center of mass, and a 
description of the object's surfaces. 

World Description: Objects, Atoms, and Transformations 

As explained in the previous section, an object's definition must include: 

(1) A mass and a center of mass for dynamical representation in constant 
acceleration fields. 

(2) A description of the object's surfaces for visual representation of the 
object at each instant of time. 

Specification of these data can be quite involved, particularly if an object is complex. 
Furthermore, objects often contain similar basic shapes: cubes, parallelopipeds, 
ellipsoids, and so forth. Accordingly, various levels for object are provided. 

The lowest level of object delmition is the atom. An atom is defined to be that 
set of data that specifies (1) and (2) as above. As each atom is defined, it is tagged 
with an atom name. Atoms form the basic building blocks from which complex objects 
may be built. 

An object is now defined to be an aggregate of transformed atoms and other objects. 
A transformation is defined to be one of the following mathematical operations 
on atom data: 

(1) RC(Cg,C{j,Cy): Recolor each surface with the color' 
C = Cg (gieen) + (blue) + Cj. (red). 

(2) C(Cg,C{j,Cj.): Color any surface that has not already been colored with 
the color C = Cg + + C^. 



Here color is additive (e.g., red + green yields yellow). 
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(3) B(b): Multiply the coefficient of optical reflectivity of each surface by b. 

(4) S(kx,ky,k2): Scale by kx in x,ky in y, and k^ in r,. 

(5) MCmj^jniyjm^): Move by m^^ in x, niy in y, and in z. 

(6) X(x,y,z,0): Rotate 6 degrees about a line parallel to the x axis and 
passing through the point (x,y,z). 

(7) Y(x,y,z,6): Rotate 0 degrees about a line parallel to the y axis and 
passing through the point (x,y,z). 

(8) Z(x,y,z,0): Rotate 0 degrees about a line parallel to the z axis and 
passing through the point (x,y,z). 

(9) W(w): Multiply the mass by w. 

This scheme of definition allows a user to build increasingly complicated objects 
out of an existing set. At no time does the user have to be concerned with what the 
basic atom data are; he merely creates new objects out of existing ones. What the user 
does, in effect, is to create a "blueprint plan" for an object; then, upon receipt of a 
display command, the objects and subobjects are constructed from the atom data in 
accordance with the blueprint plan. 

For example, suppose two atoms exist— **Cube" and **Pyramid"-and represent 
the basic shapes of Figure 11a and b. A user may then create the four-legged table of 
Figure 12 by creating a "Tabletop** out of a "Cube," a "Tableleg" out of a "Cube," 
and then "Table*' out of "Tabletop" and four "Tableleg"s. 
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Figure 12-- "Table" (original specification) 

The blueprint structure specified by the user is: 

"Table" 




Tabletop" 



"Cube" 
(Atom) 



"Pyramid" 
(Atom) 



It should be noted that an atom is always at the lowest level. Each T is a set of 
transformations and serves to link the various subobjects together. 

At this stage, the user can modify any of three objects: "Table," "Tabletop," or 
Tableleg." The highest level of definition is two. Any change that the user makes in 
one level m\\ ripple its effects into higher levels. For example, suppose the user 
redefines "Tableleg" to use "Pyramid" rather than "Cube." The data structure is now: 

"Table" 



"Tabletop" 




"Cube" 
(Atom) 



Pyramid" 
(Atom) 
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Display of "Table" now results in displaying an object with the shape of Figure 13. 
Each of the objects defining "Table" still appears, but "Tableleg" appears with its 
new definition. 




Figure /J - 'Table" ("Tableleg" modified) 



It should further be noted that this scheme is open-ended as to the level of 
definition~"Table" may now be used as easily as "Tabletop," "Tableleg," "Cube," 
and "Pyramid." For example, thfe user could create a set of several "Table"s and 
call it "Tables." 

The one thing that a user is not allowed to do is to attempt to form a recursive 
definition. Any attempt to modify "Tabletop" by adding "Table" to it would result 
in an error message since "Table" is defined, in part, by "Tabletop." 

World Description: Atom Data (Definition) 

Atoms form the data base from which all objects are constructed. The data for 
each atom contain the following information: 

(1) Mass and the location of the center of mass. (These data are used in 
the dynamics of objects.) 

(2) A description of the atom's surfaces. (These data define the visual 
properties of the atom: shape, location, and orientation in R^, color, 
optical reflectivity, and so forth.) 

Object transformations operate on and modify these data, but it is still transformed 
atom data that are decoded into a visual perspective by the image generator (except 
for color: C(Cg,C{^,Cj.) and RC(Cg,C{^,Cj.) are replacements). 

The first requirement is satisfied by encoding a quadruplet (Xg,yg,Zg,M) representing 
mass M situated at G = (Xg,yg,Zg)'^ in R^. 
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The second requirement involves encoding the boundary constraints of the type 
of [1.1] and [1.2] to allow a projection processor to project the surfaces from into 
(i.e., to project the S-dimensional world onto the 2-dimensional viewing window of 
the observer) in such a way as will then allow a visible surface processor to determine 
M^ich of the projected surfaces is actually closest to the observer (i.e., which surfaces 
are visible and which surfaces are hidden behind closer surfaces). 

This problem of projecting an object defined in 3 dimensions onto a 2-dimensional 
viewing window and representing its depth (the missing dimension) at any point on the 
viewing window is greatly simplified if the object contains only surfaces that lie on 
planes; that is, any point (x,y,z) on a surface satisfies an equation of the form 

ax + by + cz + d = 0 

where a,b,c, and d are scalar constants. Under this restriction it can be shown that 1/z 
is linearly related to (x',y'), the R^ coordinates of the projected point (see projection 
processor). Hence, a general linear algorithm can be developed for solving the hidden 
surface problem if the original surfaces are planar. 

Such planar-surface objects can easily be represented by finite sets of points, 
surfaces, and straight line edges in the x,y,z coordinate system. This is due to the fact 
that the n boundary constraints of [1.1] and [1.2] for the object reduce to a set of 
n linear equations if all of the surfaces are planar. For example, the unit volume cube 
of Figure 14 has the form 

S= { (x,y,z) I x>0,x+1>0, y>0,y+1>0, z>0,z+1>o} 
and a point on the surface must satisfy at least one of the following: 



Si: 


(x=0, 


0<y<1, 


0<z<1 




(x=1, 


0<y<1, 


0<z<1 


S3: 


(0<x<1, 


y=o, 


0<z<1 


S4: 


(0<x<1, 


y=i, 


0<z<1 




(0<x<1, 


0<y<1, 


z=0) 




(0<x<1, 


0<y<1, 


z=1) 



Each of the constraint surfaces is completely enclosed by a set of straight-line edges. 
Each straight-line edge is the locus of points that satisfy two boundary constraints. Each 
straight-line edge terminates at a vertex of the object- ia point that satisfies three or more 
of the boundary constraints. 
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Figure 14- Unit Volume Cuba 



Thus, objects with planar surfaces may be easily defined by overlaying them with a 
grid structure of points, surfaces, and straight-line edges from one object vertex point to 
another, separating one surface from another: 

Pointsf-Point data may be defined by either 

(1) Direct input of an (x,y,z) coordinate. 

(2) Specification of two of the three coordinate values and a surface. 
The third coordinate value is calculated by constraining the point 
to lie in the plane of the surface. 

Surfaces— Surfaces are defined by the input of a triplet ^1^2*^^)^ The order of 
the triplet is such that the vector N = (P^ - P2) x (PpHF^) points perpendicularly out of 
the surface. The surface is defined to be white (Cg - C^- Cj.) and has coefficient of 
optical reflectivity equal to unity. Equivalent definitions for the plane surface S of 
Figure 15 are: 

(Pc'^b'^a) (Pb'^a'^c) (^a'^c^b^ 
Edges— Edges are defined by the input of a quadruplet (Pi,P2>Sl,S|^), such 
th.iD if the edge were transversed from point P-^ to point P2, S^ is the surface on the left 
side of the edge and Sp^ is the surface on the right side of the edge. See Figure 16. 
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Figure 15 — Planar Surface S 




Non-planar analytic surfaces may be handled with a minor addition to the grid 
structure. Consider the smooth (i.e., analytic) surface Sg^ of Figure 17. The sur- 
face satisfies 

f (x) = f (x,y,z) = 0 [1.7] 
and has continuous first partial derivatives. Since f is a scalar function, the gradient 



vector is given by 



N(x) 



9f(x) 
3x 



_9f 

3x 

3f 

K 
9z 



[1.8] 
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and points perpendicularly out of the surface at each point x. Letlcg^^, and X2 be 
points on S^^ and let N(Xq), N(xj^), and N(x2) be the gradient at those points. 
Let Sfi be a planar surface passing through the points xq^^I* ^^^^2* 




Figure 17 — Smooth Surface and Flat Surface Sfg 

Since the surface is analytic in the region of interest, it may be represented by a 
Taylor*s series about the point Xq. 

f(xo+H) = f{xo) + f'(xo) (xo+h)+ -^|(xo+H)'^f''{xo) (XQ+Ti) + ^^'^^ 
= f(xo) + f'{xo)(xo+^ + o(lxo+TTl) 

where f'(xo' R{xo). 

Now let X]^ Xq + and X2 = Xq + k2h2. Clearly, as k]^,k2 ^ 0,X2,X2, xq 
and Nf(xj^),N(x2) N(xo). Hence, the locus of points on a plane passing through the 
points xq,Xj^, and X2, given by 

approaches the points on the plane 

. N(xo) = 0 

formed by the first order terms of the Taylor *s series expansion. The error is o(l x - xq I). 
Thus, the analytic surface may be arbitrarily closely a^^jximated by a set of planar surfaces. 
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All surfaces are detected by the observer as a color-brightness pair. The surface color 
is the color assigned to the object by a color transformation. The surface brightness at 
point X on the non-shiny surface is given by 



[l-lll 



B(x)« BFACT 



where: 



P 
S 

i=1 



Bgf 



lierxii iiN{x) 11/ 



JlyJI IIN(X)I 



BFACT 5= coefficient of optical reflectivity of the surface 

®AMB " brightness of ambiant light source 

P = number of point light sources 

= location of the i*^ light source in 

= brightness of the i*^ light source 

= normal to the surface at point x 

= directional vector of light from sun at <» 

= brightness of sun 

= brightness reflectivity function 
For planar surfaces N(x) = N, a constant vector, and [1.11] reduces to 



1 

Bi^ 
N(x) 

S 

BS 



[1.12] 



B(x)« BFACT 



(Bi-x)' . N 



B 



AMB 



+ S B;f 



Bcf/ . IST 



i-1 \ IIC:-xll IINII/ \ IISII IIN II 

If there are no point light sources, the brightness on a flat surface is constant and is given by 



B{x) = B 



BFACT 



BaMB + Bsf { . H 

IISII IINII 



[113] 



Now consider a non-planar analytical surface f{x) arbitrarily closely approximated by a 
set of planar surfaces. Since f(x) is analytic, f'(x) = 9f(x)/9x = N(x) is continuous. Therefore 
B(x) as given by [1.11] is continuous; in particular, B(x) is continuous across the edge 
separating one of the approximating planar surfaces from another. Since neither [112] 
nor [1.13] is continuous across these boundaries, a continuous brightness function must 
be defined over the grid structure. 
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Consider Figure 17 again. The surface corresponding to f(x) is approximated 
continuously over by a set of planar surfaces that pass through a set of points 
(xj,f(xj) = 0} known to lie on S^^. If the gradient vector N(x) = af(x)/ax is known at 
these points, then the brightness function B(x) is known at that set of points from [1.11]. 
This known set of brightness {"xj.BCxj)} allows the definition of a continuous approximating 
bri^tness function. 

Let X 6 [0,1] . Qearly, x = Xlcj + (1 - X)xj is a point that lies on the edge between 
Xj and Xj in the grid structure. Two obvious functions that could be used to approximate 
the brightness function B(x) along this edge are: 

LinearRule: B(X) = XBIXj) + (1 - X) B (xj) , X e (0 Jl [1-14] 

Logarithmic Rule: B (X) = B{xj) I Lj , Xe [0,1] [115] 

\B()?j)/ 

For maximum dynamic range for a fixed number of bits, [1.15] was chosen. 
Rule [1.15] specifies the brightness around the perimeter of each of the approximating 
planar surfaces. The brightness at any point x within a perimeter will be defined as the 
result of applying [1.15] to an edge resulting from the intersection of the approximating 
planar surface with the y = X2 plane passing through the point x = (x = xj, y = X2, 
z = X3). See Figure 18. 

Appendix A provides an example of a non-planar surface. 



x,.B(Xj) 




Figure 18 — Brightness Within Planar Surface Perimeter 



FinaS Atom Data 

Once the grid structure over an atom has been completed, the atom is preprocessed 
and put into the atom library in a compact form. This preprocessing involves: 

Joints: All point values are calculated and points are stored as an x = {x,y,z)^ 
96-bit floating point ta-iplet. 

Surfaces : All surface normals are calculated to insure that the (Pi,P2'^3) 
triplet does, indeed, define a surface normal N = 

Brightness Normals: Brightness normals are calculated for vertex points on 
analytic, non-planar surfaces. These normals are held to 10-bit precision per component. 

Edges : The general form of an edge becomes 

E = (Pi.P2'P3'P4'Sl'Sr,Bi,B233,B4) 

where 





= label of end point Pj^ 


P2 


= label of end point P2 


P3 


= label of point on left surface 


P4 


= label of point on right surface 


Sl 


= label of left surface going from Pj^ to P2 




= label of ri^t surface going from Pj^ to P2 



81,823334 appear only if Sl or Sj^ is a non-planar surface. If Sl and Sj^ are approxi- 
mating planes for the same analytic surface, B-^ is the label of the bri^itness normal at 
Pj, and 82 is the bri^tness normal at Otherwise, 8-j^ and 82 refer to bri^itness 
normals for Sl at Pj^ and P2, and 83 and 84 refer to brightness normals for Sj^ at 
?l and P2. 

The atom data remembered in the atom library and used in decoding the image are: 

(1) Point data 

(2) 8rightness normal data 

(3) Edge data 
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Section 2 
PROJECTION PROCESSOR 



INTRODUCTION 

This section provides a detailed description of the projection processor of the 
CHARGE image generator. The projection processor accepts as input a description of 
a 3-dimensional space from an observer within this space. The nature of this description 
is detailed in Section 1, Worid Description and Viewing Parameters (page 20). The 
projection processor performs two functions: 

(1) It reduces the total set of 3-dimensional world description data into 
a subset that is potentially visible to the observer. . 

(2) It transforms this set of potentially visible 3-dimensional data into an 
equivalent set of encoded 2-dimensional potentially visible data. 

The latter data set serves as input to the high-speed linear visible surface processor that 
forms the video image and then encodes it for display on a CHARGE display terminal. 

Input Data 

Input data to the projection processor contains the following information: 

(1) Specification of the observer in R^ 

(2) Specification of the light sources in R^ 

(3) • Mappings of atoms into R^ 

Observer 

The observer is completely specified by the set of parameters (XY^r,Y^,Z;y,X!j;, 
Y'e>Z'e,Xy,Yy,Zy,W'j^,W'y,MAG,Bo) which define an eye, a viewing window, and a 
point of observation. These parameters are pictorially represented in Figure 19 
V = XyjYyjZy)'^ defines the point of observation with respect to "O = (0,0,0)'^, the 
origin of R^, This is the point toward which the observer directs his vision. The observer 
sees mto R^ only through a viewing window. This viewing window has a center located 
at W = (X^,Y^,Z^)'^ with respect to the origin of R^, and lies in a plane orthogonal to 
the vector V - W. 
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W is the origin of a second rectangular axis coordinate sj . with axes x',y',z' and 
unit directional vectors TJ', and k' defined so that: 

7x7«o 7x? = -r k'xT«7 

r X r « r r X r = 5 k' x y « -t^ 
r X TJ' « -T r X F = r f x ij' » is 

Ti (V' - W) 



k 



II V' - wll 

7 .T = Ty 



Mth respect to this primed coordinate system, the window size is W'^ horizontally 
by Wy vertically, and the y and y' axes lie in a plane. The observer's eye is located at 
■ (X'E,Y'g,Z'£j)T with respect to the primed axes. If X'j; = Y'e = 0, the eye looks 
at V - (Xy,Yy,Zy)T directly through W = (X^,Y^,Z^)T the center of the viewing 
window. The observer looks through the window with a telescope with magnification 
factor MAG > 0 and coefficient of optical transmission B^e [0,1] . The relationships 
between the observer, the viewing window, and the point of observation are pictured 
in Figure 19. 

Light Sources 

The types of light sources are provided for: a planar light source and point li^it 
sources. The planar light source is specified by defining a sun vector and brightness 
level fraction 

L$ = (Xs'Bs) [2.1] 

where Xg » (^tg.yg.Zg)'^ specifies the durection of the parallel light rays through R^, and 

€ [0,1] specifies the brightness as a fraction of Bj^^^^, the maximum allowed 

brightness. This sunlight source illuminates a "non-shiny" surface with gradient vector 
— 9f(x) , 
N(x) = _ with brightness 
9x 

xJR(x) ^ 

117,11 IIN(x)li; ^^^MAX [2.2] 

where f(a) is a reflectivity function. All surfaces with the same gradient vector will be 
illuminated with the same brightness. 
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Figure 75- Observer in 

A point light source is specified by defining a location vector, a brightness level 
fraction, and a brightness fall-off fraction. 

L- (xi.Bi.r;) [2.3] 
where xj = (xi,yi,zi)'^ specifies the location of the i*^ point light source of brightness 
^i%AX ^ f ^ ''i- P°'"* source emits light radially from xj with 
brightness diminishing inversely with the square of the normalized radical distance 
(f/r ) from xj 



lB(r/ro.)l = 



Bi^MAX 



[2.4] 
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This point li^t source illuminates a surface at point x with gradient vector (N(x) - 
3f(x)/3x with brightness 

\llx-i-xll liNBlll/ /||i;,-;|| y 



as shown in Figure 20. The world description contains only one sunlight source, but may 
contain an arbitrary number of point light sources. 




Atoms 

Atoms, the basic building blocks of objects in the world, are mapped into their final 

q 

orientations in R*^ by composite transformations that link object definitions from their 
highest level down to the atom level Transformations that map the atom's point data 
in are implemented via 4x4 matrix operators. These are: 



0 0 

1 0 
0 1 
0 0 



m. 



[2.6] 
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S{kj^, ky, kjj) 



0 
0 

Lo 



0 

ky 

0 
0 



0 
0 
k. 



0 
0 
0 
1 



[2.71 



Y{x,y,z,OJ 
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sinOy 


0 


COS(?y 


0 


0 


0 


0 


1 _ 


cosd^ 
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cosO^ 
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1 
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1 « 



M{-x,-y,-z) 



[2.81 



M{-x,-y,-z) f2.9] 



M {-x,-y,~z) 



[2.101 



The composite transformation [T*l that maps an atom in an object of level of defini- 
tion m into its position in is given by 



iri = lTJlT^_ll ...IT3IIT2IIT1I 



[2.111 



where [Tj^] is the transformation (which may itself be composite) which defines the 
atom*s orientation in the k level of object definition in terms of its orientation in the 
(k- 1)* level of object definition. 

[T*l maps a point = (Xo>yo»2o>l)'^ original atom definition into its final 

position Pf = (Xf,yf,Zf,l)'^ via the matrix product 

Pf = [ripQ [2.121 

The atom definition is such that its final form in is strictly a function of its points 
in R^. Hence, [T*l : Atom (original) Atom (final orientation in R^). 
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Processing Steps 

The projection processor projects the 3-dimensional world description 
^Observer. Sunlight Source. (Pt. Light Sourcosj , {{[Tj] , Atomj)| J 
into an encoded 2-dimensionaI description via the processing steps of Figure 21. 



Add R 2 Atom 

Data To 
Output Data Set 



YES 




Project Atom 
Onto r2 Viewing 
Window 




NO 




YES 



Y Order the 
Output Data 
Set 



STOP 



NO 



Figure 21 — Projection Processor Processing Steps 
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The amount of processing is roughly equal to the number of atoms in the world 
description times the average atom complexity. The output data set that will be passed 
to the visible surface processor may contain data from all of the atoms or, at the other 
extreme, it may be empty. 

Initial Processing 

The job of the projection processor is to encode the world description onto the 
viewmg window in accordance with the observer's parameters. This is most easily 
accomplished if (a) the observer's eye is located at the origin 0 of R^; (b) the observer's 
viewing window is located in a constant z plane with the x' axis parallel to the x axis, 
the y' axis parallel to the y axis, and the z axis parallel to the z axis; and (c) the 
observer's apparent magnification factor is unity. Let [T^] denote the matrix transforma- 
tion required to do this. If [T^] is applied to the observer, the observer's viewing window, 
each of the light sources, and each of the atoms in R^, then the relative positioning of 
this transformed input data set is equivalent to that of the original inpat data set. 

Transformation [T^] may be computed in two stages from the observer's parameters. 
The matrix transformation 



where 



IT'wIMTmAgI [TxHTyllT^il 

JTmll : M(-X^,-Y^^,-Z^) 
(Tyl : Y(O,O,O,0y) 0y = tan-1 



(Xy Xyy) 



^l (Zv - Z^)2 



[Txl : X(0, 0,0,0^) = tan"l 

^l(Xv- X^)^ + (Zv - Z^)2 

[T|V|agI- S(MAG, mag, MAG) 



moves the observer and his viewing window from the orientation of Figure 19 into that 
of Figure 22. The primed and unprimed axes now are coincident with the center of the 
window now lying at 0, the origin of R^. The viewing window now lies in the z = 0 
plane and the eye is located at E = (Xe,Ye,Ze)'^ where = MAG-X'e, = MAG-Y e, 
and Zg = MAG-Z'ij. If the viewing window is left the same W'^^ by W y size, the desired 
magnification is accomplished. 
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il j 

w=o 




Figure 22- World After [T' ] Transformation 



The final desired magnification may be obtained by applying a move transformation 
of M (- Xg,- Yg,- Zg). This transformation applies to the window coordinates as well 
this time and results in the observer-viewing window orientation of Figure 23. 
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MAG • Y'e 












X4 








E=0 


>^ 


mag-x'e ^ 




^ 



2 A 

Figure 23 - World After [T,„] Transformation 
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And so 



[Twl = [Tm2] [TmaI 1 [Txl [Tyl [Tmi 1 [2.13] 

where 

(TmiI : M(-Xw, -Yw, -ZJ 

(Ty] : Y(O.O,O,0y) Oy = tan"! 



\f(Z^- Z 



(Txl : X(0,0,0,0,) 0, = tan'l "'Yv Yw) 



\|(Xv- X^v)'^+ (Zv- Zj2 

ITmagI : S(MAG,MAG,MAG) 

[T|vi2] : M(-MAG-X'e, -MAG-Y'e, -MAG-Z'e) 

The remaining window now has comer coordinates: 
W'x W' 

Xi = (-— -MAG-X'e, MAG- Y'e, -MAG-Z'e) = (X'a/||N. Y'maX. ^vw) 

_ w' w' 

X2 = (+-2^- MAG-X'e,+^-MAG- Y'e,- MAG-Z'e) = (X'maX- Y'maX. ^vw) 

f2.14] 

X3 = (- W;,-MAG -X'e, - W'y - MAG- Y'e, - MAG-Z'e) = (X'min- Y'^iN, Z^w) 
2 ~ 

- w'x w'„ 

X4 = |+-^-MAG- X'E,-yl-MAG- Y'e,- MAG-Z"e) = (X'min. Y'min, Z^) 
and the initial processing steps are: 



/' \y' -7' 







(Xv, 


[Twl / 






— Ill 


[Twl 


Xsll 



,0. . - [Twl / [TwlXs \ 

(2) L3 = ,X3,B3) (^-^^pj:^ ,83) = (S„,B3) [2.14a] 

(3) Lj = (Xj.Bi.ro.) !H ([TjXj.Bj.MAG-ri) = (^^i. Bj, ro'. ) 



As was mentioned, [%] must also be applied to ([TjJ.ATOMj) pairs to maintain the 
original relative positioning of the viewer, light sources, and objects. Functionally 

[Til [T^^] 
Atom in atom library »■ Atom in world y Atom fn transformed world 

This latter operation of transforming l[Ti],ATOMi) pairs into ([T^] [TiJ.ATOMj) pairs 
will be done in the atom projection loop. 
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Atom Projection Loop 

Object data entering the projection processor consists of ([T^],ATOMi) pairs where 
[TJ is a matrix transformation that maps the points of atom ATOMj into (original) 
position. Also available is matrix transformation [T^] that maps the original into 
an equivalent 3-dimensional space in which the center of the viewing window is at the 
origin of this new space. Atom data consists of (a) a finite set of points defining the 
atom*s vertices in its standard (library) orientations, (b) a finite set of brightness normals 
representing the gradient vector 3fj(x)/0x of the atom*s surfaces in the atom*s standard 
(library) orientation, and (c) a finite set of straight line atom edges which link together 
the atom*s surface constraints and brightness functions, the edge definition is invariant 
to the atom's orientation. 

The first processing step is to transform the atom into its final orientation. Points 
are transformed via [2.15] . 

Pf = [TJ [Til Pq [2.15] 

where [T^] [TJ is the composite transformation mapping the point P^'^ = (^o'^o'^o'^^ 
into its final position Pj'^ = (Xf,yf,Zf,l). The set of surface normals {9fi(x)/9x| are 
transformed via [2.16] . 

Nf = AdjT ([TJ [Til} No [2.16] 

where is an initial gradient vector Nq'^ ^ (nx,ny,n2). [2.16] may be easily derived 
as follows: Let N be orthogonal to surface S and let x^ be a point on S. 
Then jvjT . (jj- ^ ^ q ^-^g^ x ^ Xq 

Now let [T] be a composite matrix transformation mapping into R^. Let [T*] be 
another matrix transformation such that 

([T*]N)T. ([T] (X- Xq)) - 0 VxeS, x^Xq [2.17] 

That is, [T*] maps N in such a way that the transformed normal vector [T*] N is still 
orthogonal to the transformed surface. Now since ([T*] N)^ = n'^[T*]T, [2.17] may 
be rewritten as 

n"^[T*] "^[T] (X - Xq) - 0 V X S, X 9^ Xq [2.18] 

Now [2.18] is obviously true whenever 

Un'^lT] - k[l] [2.19] 
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where k is any non-zero scalar. Now, if [T] is of the form 



[T] = 



[A] 



3x3 



[01 



3x1 



"3x1 



1x1 



it has an inverse only if [AJg^^g is nonsingular. Brief examination of [2.6] • [2.10] 
reveals that all simple matrix transformations are of this form, providing that the scaling 
transformation S(kj^,ky,k2) has all components > 0. Furthermore, composite trans- 
formations are of the type 



(Tl 



n 






n 


[Tjl = [Til 


•IT2I . . 


i=1 






"n 




n i-1 


2 


tAil3x3 ; 


n (n 


i=1 




i=i i=i 




[01 1x3 1 





[lllxl 



where [Tj] is a simple matrix transformation and so [T"^] is guaranteed to exist. 
Thus, [2.19] is equivalent to 



ITI 



Adj[T] 



or, equivalently 



[2.20] 
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which proves [2.16] . 

It should be noted that the brightness functions [2.2] and [2.5] involve normalized 
gradient vectors N/||N||. Unfortunately, [T*] as given by [2.20] is not, in general, a 
unitary transformation (due to possible unequal scalings k^^.ky, and k2) and so the set of 
gradient vectors |Nj(x) = 9fj(x)/9x| must be normalized after application of [2.15] . 

Edge data containing only point labels, surface labels, and gradient vector labels 
are invariant under any composite transformation. 

Projection of into R^ : Consider Figure 24. The point P = (x,y,z)'^ is in R^. 
The viewing window is located in the constant 2 = Z^^ plane. The eye is located at the 
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origin O! The projection point P' = (x',y', Z^^)'^ of the point P onto the viewing window 
is desired. 



P=|x,y,z) 




Figure 24 - Projection of Point P Onto the Viewing Window 



Now, since P' = aP, where a is scalar, 



X 

V 

z' = Z, 



VW 



ax 
ay 



and from the last, a = Zyy^jz, That is, a is the ratio of the viewing window depth to the 
total depth of the point. With a so defined, the first two equations reduce to 



-vw 



-vw 



y = 



[2.21] 
[2.22] 



Both [2.21] and [2.22] have singularities at z = 0. That is, points on a plane that are 
parallel to the viewing window and that also contain the eye have no valid projection. 
Also, points with z < 0 have no valid projection since they lie "behind" the eye. 

Note too from [2.21] and [2.22] that the ratios x'/x and y'/y are greater than 
unity if 0 < z < Z^^. That is, an object wil) be magnified if it lies between the eye and the 
viewing window. Because of this magnific.tion and the singularities of [2.21] and [2.22] , 
a restriction will be made that objects lie on the "far** side of the viewing window, that is, 
objects will not be allowed to penetrate tht viewing window plane and still be seen. If 
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a surface does penetrate the viewing window, only that part of it which lies on the far side 
of the viewing window will be displayed. 

Now consider the projection of a straight-line edge from into as in Figure 25, 
Since Pi,P2. and E define a plane containing Pj^ - Pg which intersects the plane of the 
viewing window (z = Z^^) in a straight line, [2,21] and [2,22] map straight lines into 
straight lines. 




(PrP2) 



Figure 25 - Projection of a Straight-Line Edge 



dx 



The slope m = — of a projected straight-line edge may be found as a function of the 
original R^ point data, m is given by 



m 



X 9 - X 



y2" y 1 



where (x'i,y'i,Zy^y) and 2''^v/^ correspond to projected points P'j^ and P'2- 

Application of [2.21] and [2.22] yields 



-vw 
^2 



1X2 



m 



-vw 



{——) 72" ( )yi 

^2 ^1 



Z1X2Z2'«1 



m ' 



2172 ziVi 

Note that m = Adx'/dy' is the inverse of its normal definition and that here 

Ax' - mAy' 



[2.23] 



[2.24] 
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and, hsWe, that the locus of points on a straight-line edge that has been projected onto 
the viewing window has the form 

X (7 2) = x' (/-,) + m(y '2-/1) [2.25] 

where (x'(y'i),y'i) and (x (y 2),y 2)^ Note also that m 0/0 if both zj^,Z2 0. 
That is, there is no valid projection of edges that lie in the eye plane parallel to the 
viewing window plane (i.e., the z = 0 plane). 

Now consider the projection of a planar surface from into R^. Since any 
surface on an atom in R^ is a planar surface bounded by straight-line edges, the area in 
R^ corresponding to that surface is that area bounded by the projected edges. See 
Figure 26. 




Figure 26- Projected Atom Planar Surface S 

Now consider whether or n^i surface S will be visible after it has been projected 
into R^. Three conditions are obvious: 

(1) The surface must face toward the eye and not away from it. Let N 
be the gradient vector of planar surface S. Then this condition may 
be tested at any point P on S by 

' (P-E) = n"^ • P<0 [2.26] 

(2) Some part of the projected surface S' must lie within the region 
I (^'^y') min ^ ^' ^ ^'max»y'min <y < y'max 1 where 
^'min»^'max»y'min y'max define the rectangular boundaries of 
the viewing \vindow in the z = Z^^ plane and may be found 
from [2.14]. 
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(3) Some part of the surface in R'^ must lie on the "far'' side of the 
viewing window. Thus restriction was placed after developing the 
projection algorithms [2.21] and 12.22]. If all the edges bounding 
surface S in He on the "near" (i.e., z> Z^^) side of the viewing 
window, then S' will not appear. 
A surface S that meets these three conditions is potentially visible over the 
region { (x',y') I (x',y') e S'} on the viewing window. Whether or not it will be visible 
at any position (x',y') eS' depends upon whether or not S is the "nearest" potentially 
visible surface encountered by an eye sighting from E = 0 through that point 
(x',y', Z^^)'^ on the vie\ving window. This distance between a surface and the viewng 
eye will now be related to the projection of the edges that define the "right" boundaries 
of S\ 

First, consider the relationship in R^ between projected edges and surfaces. 
The original atom edge definition (Pi,P2»P3»P4,Sl,Si^,Bi,B233,B4) is such that if one 
were to travel from end point P-^ to end point he would encounter surface Sl on 
the left side of the edge and surface Sj^ on the right side of the edge. If the edge is 
projected by [2.21] and [2.22] into R^, there are four possible outcomes: 

(1) Both Sl and Sj^ face away from the eye; that is, both surfaces 
fail to meet Condition (1) for a surface to be potentially visible. This case is that of 
an edge on the far side of an object.. The edge is ignored. 

(2) The edge projection in r2 is horizontal. In this case, one projected 
surface is above the projected edge and the other projected surface is below the projected 
edge. Since neither projected surface has any left or right orientation to the edge, the 
edge is ignored. 

(3) -At least one of the surfaces meets condition (1) and P'2> the pro- 
jection of point P2, lies "above" P';^, the projection of point P^. In this case, the 
original left-right convention remains valid. If either of the surfaces does fail to meet 
Condition (1), its surface label S is replaced by S which will flag to the visible surface 
processor that only the other surface is potentially visible. Figure 27 illustrates this case. 

(4) At least one of the projected surfaces meets Condition (1) and P'2 
the projection of end point Pg, lies "below" P'^. In this case, thp left-right convention 
is no longer valid. This may be corrected by interchanging all of the labels as shown in 
Figure 28. This case then becomes the same as the third outcome. 

The relationship between the depth z of a surface in R^ and one of its "right" 
boundary edges will now be developed. Consider Figure 29. P2 =^ (x2,y2»Z2)'^ been 
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f/5t/re2S- Projected Edge With P'l Above P'2 




determined to produce a projected point P'2 which lies ''above*' P'j^. P = (x,y,z)^ is any 
point on planar surface S. The eye, E = {x^.y^^z^)^, is located at the origin 0. The 
viewing window lies in the z = Z^^ plane and F = (x',y', Z^^)'^ is the projection of P 
onto the viewing window. N = (n^^nyjn^)"^ is the gradient vector for surface S, 

Now, since both P2, and P are on surface S, the vector P- P2 lies on S and, 
hence, is orthogonal to gradient vector N. Thus, 

(P-P2)T.N= 0 [2.271 
which may be expressed as the scalar equation 

nx(x - X2) + ny(y - + n^lz - Z2) = 0 [2.28] 

Also, [2.21] and [2.2] may be solved for the x and y coordinates of the original point 
P in terms of the z component of P and the projection point P'. 

X = (z/Zy^lx' 
y = (z/Zv^)y' 

substituting into [2.28] and multiplying Zy^ yields 

nxlz • x' - Z^w • X2) + ny (z • y' " Z^^ * 72) + n^^Z^^ • z " Zyw * = 0 
zfRT • P') - Zv^InT . P2) = 0 

or simply 

d = — = —— [2.29] 
2 N ' . P2 

For surfaces lymg on the "far** side of the viewing window, z ranges from Zy^y to - 
and maps d from +1 to 0"**. Thus, ze[Zy^r^) implies de[l,0). Thus, the problem of 
determining the nearest surface seen by an eye that is located at the origin and sights 
through the point P' = (x',y',z' = Z^^) is that of determining the surface with the 
maximum d. ' [2.29] can be written as 



d = 



n^x' + n y' + n^Z 



z vw 



K 



d = 



"z^vw "x 



+ — X + — y 
K K 



[2.30] 
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where [2.30] is a function of the position (x',y') on the viewing window. 



d. = JJ^ + Ji x'^ y'o [2.311 
° K K 2 K ^2 



Thus, if x' = x'2 + Ax' and y' = y'2 + Ay', then 

d = -iZl! + JL (x' + Ax') + (y'^ + Ay' ) 
K K 2 ' K 2 y 

= do + Ax' + Ay' 
K K 

d = do + dxAx' + dyAy' [2.32] 

where d^^ = ^xl^ dy = ny/K. Thus, [2.32] gives the depth function d for a displace- 
ment of (Ax , Ay ) in R . Consider now a displacement in R"^ along the projected edge. 
If the slope of the edge is m, a movement of Ay' vertically yields a horizontal movement 
Ax' = mAy' and the depth function is given by 

d = do + dx(mAy') + dy (Ay') 

= do + (dxm+dy)Ay' 
d = do + djAy' [2.33] 

where dg = md^^ + dy. Thus an expression for the depth function of a projected planar 
surface where the displacement in R^ is the composite of Ay' vertically along the projected 
edge, and then Ax' horizontally from the projected edge is given by [2.34] and is 
illustrated by Figure 30. 



+ dxx' + dgAy' 



[2.34] 



where: 
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''2' ,nx, , ( "^"x ^ "y , [2.35] 
K K K 



n'^ = (n3^,ny,n2) = Gradient vector to surface S 

= {^2>y2'^2^ ~ "Top'' end point of edge 
P'2'^ f^'2'y'2'^vw^ ~ Projected top and point of edge 
K = • = nxX2 + nyyg + n^Zg 
m - Slope of projected edge 
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Ay' = (y' - y'2) = Vertical displacement in R"^ from P'2 
Ax' = ^x' - (x'2 + niAy'j^ = Horizontal displacement in from edge. 
Note that in Figure 30 both Ay' and Ax' are negative. 




X 



Figure 30- Depth at Arbitrary Point P 
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The final function defined over R^ which is to be projected onto the R^ 
viewing window is the brightness function B(x) given by [2.36] (derived from [1.11] 
and [2.14a] .)• 

[2.36] 

f /(xrx)^.N(x) 




llxi-xll llN(x)ll 



where: 



B^: Coefficient of optical transmission for observer; Bq€ [0,1] . 

^max' Maximum allowed brightness of any one source; Bj^^^ > 0. 

Bp: Coefficient of optical reflectivity of surface; Bpe[0,l]. 

f(a): Reflectivity function. 
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Bg: Fraction of B^^^ that serves as brightness of sun light source; 

Bs€[0,ll. 
S: Direction of tlie sun's light rays. 
N(x): Gradient vector to surface at point x. 
P: Number of point light sources. 

Bj: Function of Bj^^^c ^^^^ source's maximum brightness; 

BiG[0,l]. 

5C|: Location in R of the i*'^ light source. 

roj: A constant such that 11 xj - xH /rg. is a brightness fall-off ratio for 
the i*'^ light source; rQ. > 0. 

There are several possible brightness definitions that an edge may have: 

(1) Both the left and right surfaces are "flat." For this case, Nl may 

serve for both points ?i and ?2 ^^^^ surface and Nj^ may serve for both and 

P2 for the right surface. 

SPECIAL CASE: If (1) holds and if P = 0. then the brightness of each 
surface is a constant independent of x (see [1.13]. These surfaces will 
bo called flat-shaded surfaces, 

(2) The edge is to be smoothly shaded across (that is, B(x) is continuous 

across the edge) and N = (Nj^ + %)/llN2 + %ll defines the gradient vector at any point 

X on the edge. In this case there is one brightness at ?i and a second brightness at 

SPECIAL CASE: If (2) holds and if P = 0. then the brightness along 
the edge is constant. 

(3) The edge is to be smoothly shaded across with Nj^ defining the bright- 
ness normal at Pj^ and N2 defining the brightness normal at P2. lliis is the typical form 
of an edge on a planar surface facet that is used to approximate a curved surface over 
some region. 

(4) The edge is either flat-to-curved, curved-to-flat, or curved-to-curved. 
These edges are typically the result of the intersection of two atom surface constraints, 
for example, the junction of the flat ends and the curved sides of a right circular 
cylinder. The brightness function is discontinuous over these edges. Each of these 
edges generates two projected edges, one for the left-side surface using Bl and B2, and 
one for the right-side surface using B3 and B4. 

Thus, the brightness of an atom edge is typically specified at its end points by 
application of [2.36] to the brightness normals at the .'Appropriate end points. Since 
Case (4) generates two projected edges, the general case may be taken to be that of 
Figure 31 without loss of generality. 
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Figure 31 - Projection of an Edge With Two Brightness Normals 



The continuous brightness function along the edge had been defined as the logarith 
function (see [1.15]). 
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where X e[0,l] represents the fractional distance of x" from Pg to P^; that is, 
P = X(Pi.) + (1-X)P2. However, P' + (1-X)P'2 and, so, all smooth shading will 

be done in the perspective domain approximating the shading in R^. 
Thus, since in R^^an edge has the form 

x'(y') = x'ly'j) + m (y' - y'2) 
which is linear in y', the approximating brightness can be written as 

y ^ y 2 

■B(yV\''^^ 



B(y') = BIy'j) 



B(y', 



where 



B(y'2) = BIP'j) = BIPj) 
B(y',) - B(P',) = B(P,) 



[2.37] 
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For computational purposes, the form of [2-37] may be altered. Taking the logarithm 
of both sides of [2.37] 

(y' y' ) 

'09a (B(y')) = ■oga(B(y'2))+ "y^jT^ K (bIvV) ' (^^^2^)] 
or, more simply 

b = + (y- y'jlg [2.381 

where 

= logg (B(y'20 = 2n B(y'2) / 2n a [2.39] 



[loga(B(y'i)- logg (Blyg))] [gn B(y'^)- gn Bly'gj] 

(y'l - yiz) " (y'l - y'2)2n a 



[2.40] 



Thus, the brightness function along an edge may be put into a form identical to that of 
[2.33] , the depth function d along the edge. 

The video output at any point P' = (x', y', z^^)'^ on projected surface S is given 
by the vector 

V(P') = B(P')C'(S) [2.41] 
where C'(S) = (c'g, c'^, c'^)'^ defines the transformed color of surface S. Then 

loga(V(P')) = logg (B(P') C'(S)) [2.42] 

or, more simply 

7(P') = b(P') + c(S) [2.43] 

where V and care logarithmic vector functions and b(P') = (b(P'), b(F'), b(P')) 
Equation [2.43] has the scalar form 

green (P') = b(P') + Cg 

blue (P') = b(P') + Cb 

red (P') = b(P') + Cr [2.44] 

where Cg, c'b, and Cj are logarithmic measures of c'g, c^, and c'^. Of course, the terminal 
must perform the antilog functions 

V(P') = log;l(v(P')) [2.45] 
or, for each of the three color components, 

b(P')+Ci 

Vj(P) = a ' i= 1,2,3 [2.46] 
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This is most easily done by giving the three video amplifiers an exponential (i.e., anti- 
logarithmic) response. 

The value of logarithmic base "a" is calculated in Appendix B. 

Window Clipping: X and Y 

The three requirements for a plane surface S in to be potentially visible in 
r2 are: 

(1) The surface must face toward the eye. 

(2) Some part of the projected surface must lie in the region ((x',y') I x^^^ < 

^' ^ ^'max' y min ^ y' ^ y'maxi ^^here x'j^j^, x'^^, y'^^^^ and y'^^x 
define the rectangular viewing window in the z = Zy^ plane, 

(3) Some part of the surface S in must lie on the far side of the 
viewing window. 

Condition 1 is satisfied when edge parameters are projected onto the plane of the 
viewing window. That procedure is described in Section 1, Atom Projection Loop (pg. 45). 
This section describes implementation of Conditions 2 and 3. 

The first step involved in the clipping process is to reference the projected edge 
parameters to a more convenient set of axes. Consider Figure 32, The original viewing 
window is shown relative to the original x' and y' axes and z = Z^^. A new set of 
axes X and Y is drawn such that X increases positively to the left of x'j^^^^* and Y' 
increases positively downward from y'max- 

Also corresponding to the viewing window is a resolution grid structure of Rx 
elements in x by Ry elements in y. For stan^iard television Rx=640 and Ry=480 yielding 
the standard aspect ratio Rx/Ry^4/3, For correct (non-distorted) imaging, the observer's 
viewing window should be specified such that Wx/Wy==4/3, For the CHARGE terminal 
Rx=1600 and Ry= 1200 yielding Rx/Ry=4/3 as its aspect ratio if Wx/Wy=4/3. For 
non-standard video display formats, the grid resolution should reflect the Wx/Wy ratio. 
For example, suppose only the top-half of the television screen is to be used for a 
particular image. Then W'j^/W'y=4/1,5. For the resulting image to be unwarped, 
R'x/R'y=4/1.5 or R x=1600 (for full horizontal image on the CHARGE terminal) 
and R'y=600, The resulting display would have projected objects of the same size and 
shape as the original display, but the top and bottom quarters of the original display 
would no longer be visible. 

Figure 33 demonstrates one po.ssible implementation of Wx/Wy "4/1 .5. The unused 
scanlines (l/2(1200)-600) were defined to be below the used scanlines, these unused 
scanlines may then be used to display a second image concurrently with the first. 
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Figure 32- Translated Axes X and Y With Corresponding Viewing Window 
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Figure 33- Effect of Wx/Wy^4/3 (Rx/Ry - W^/Wy 
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The convention that will be used throughout the remainder of this section is that 
the viewing window is units (meters) wide by Wy units high and that the image 
is to be displayed over a grid structure of resolution points horizontally by Ry 
scanlines vertically. For the CHARGE termmal, Rx < 1600 and Ry < 1200. With this 
convention, the line x' =^'max "^^PS into the line X = 0; the line x ^x'^^i^ "^^Ps i^^to 
the line X = Rx; the line y' = y'max "^^Ps into the line Y = 0; and the line y' = y min 
maps into the line Y = Ry. 

The general rule for mapping points is 



X = 



max ^ ' 



Xo + ax' 



max 



- y') = + by' 



[2.47] 
[2.48] 



Slopes of projected edges are now given by 

_Ax_ _ (Xp+ax'^) - (Xp+ax'g) 

Ay " (Yjj+by'i) - (Yjj+by'2) 

a(x'i-x'2) a 



M = 



b(y'i-y'2) b 

"x \^y 
y X 



m 



[2.49] 



If (Rx/Ry) - (Wx/Wy), as is always the case for an undistorted picture (partitioning of 
the window only), then 



M 



(W^/Wy) 



= m 



[2.50] 



The incremental depth function along a horizontal line is given by 



D„ = 



D(Xi) 


- D(X2) 


Xl 




d(X'i 


) - d(X'2) 


Xl 


-X2 


d(X'i 


) - d(X'2) 


a(X'^ 


- x'2) 



7 



[2.51] 
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The incremental depth function along a projected edge is 



D(Yi)- 


D(Y2) 




^2 


d(y'i)- 




1 




d(y'i)- 


d(y'2) 


b(y'i 


-y'2) 











[2.52] 

Obviously, B(X,Y) = B{x', y'). The gradient of brightness along an edge is now given by 

[loggBIP'-,)- loQgBIP'j)! 

G = Y 9 12-53] 

All of the edge parameters have nov/ been converted to the new axis system. 

Example: Let Wx = 1/3 meter, Wy = 1/4 meter and let Rx = 1600 and 

Ry = 1200 (full screen display). Then Wx/VVy = 4/3 = Rx/Ry (image is undistorted). Also, 
a = -Rx/Wx = -4800 resolution pts/meter 
b = -Ry/Wy = -4800 scanlines/meter 
m = k meters/meter = > M = k resolution pts/scanline 

dx units/meter Dx = ^gQQ units/resolution pt 
■ds 

dg units/meter Dg = ^gQQ units/scanline 

- G 

G units/meter g = "ToKfT units/scanline 

Constants that are valid for the entire image and which may be evaluated during initial 
processing are: 

W 

'^'max -f - MAG • X'^ 



W 

y max ^ -f - f^AG ■ Y'e [2.54a] 



a = -Rx/W>, -> a"'' = 1/a 

-> c = a • b" ^ 
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b = -Ry/Wy -> b"'' = 1/b 
^0 ° ^^^'max 



max 



Initial processing of atom edge parameters yields: 



[2.54b] 



K = • P2 = • X2 + Oy • 72 + • Z2 
Pi - (Xi,yi,Zi) -^^^ x^, Yv^vwj 



P2 (X2'y2'22' I 7^ ^2' 7« y2'^vw 



[2.551 



m = 



! 2 X 1 ^1 * ^2 ~ ^2 * ^1 



y 2 - y 1 • y2 - Z2 • y. 



/ S"^ • N(Xi)\ P B, 



a s 



Jlxj-XjII llN(Xj, 



i=1,2 



where K also serves as test [2.26] as to whether or not the surface faces toward the eye. 
Final processing of the edge parameters yields: 



by'i + 



M = cm 



i = 1,2 
i = 1.2 



-1 "X -1 J 



(mn + n^) 
D3 = b-1 .-JL_VL=b-l -(13 

b(Yi)= log3(B(P'i)) 

g = b-llblV,)- bIYj)] /(y'l - y'2) 



i= 1,2 



[2.56] 
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Now consider the transformed projection of an edge onto the z=Z^ plane. The 
viewing window R is {(X,Y) I 0 < X < R^, 0 < Y < Ry | . The goal of the X-Y clipping 
unit is to restrict the set of transformed projected edges to. a subset that might generate 
projected surfaces {S'} such that S'^GR. There are several possible cases vertically: 

(1) If Yj^< 0, the entire projected edge lies above the viewing window 
(since Yj^). Since an edge serves as a boundary between one surface on its left side 
and another surface on its right side, the edge may be deleted from the projected atom 
edge set. 

(2) Similarly, if Y2> Ry, the edge lies entirely below the viewing window and 
the edge may be deleted from the projected atom edge set. 

(3) If an edge extends below Y = Ry, it may be arbitrarily defined to end 

at Y = Ry. 

(4) If an edge starts above the viewing window but extends into the region 
1 < Y < Ry, it may be arbitrarily started at Y = 1. This means 

Y 2 = ' 

X'2 - X2 + M(l-Y2) 
Do = + D3{l-Y2) 
^top = btop + g(l-Y2) 

Figures 34 and 35 demonstrate how edges are clipped vertically into the proper Y range. 

Edges that survive the vertical clipping operation are clipped into the proper hori- 
zontal X range. It should be remembered that edges serve three functions, (a) they form 
left-right partitions between projected atom surfaces, (b) they define the depth functions 
and the color of the left-side projected surface, and (c) they define the brightness along 
the perimeters of projected atom surfaces. Thus, if the set of projected atom edges that 
cross a Scanline Y^ were to be scanned in order of increasing X, each edge encountered 
would serve to (a) "end'' the surface on its right side, (b) "start'' a surface on its left 
side, and (c) define a known brightness to be used by the logarithmic brightness 
function [1.15] within the surface perimeters. These are operations that the visible 
surface processor will have to perform. 

The first operation, X-ordering a set of projected edges, is impossible to accomplish 
via one-pass techniques if the range of possible X values is an mfinite set. Furthermore, 
the only surfaces that will be seen lie in the interval Xc[0,Rx] . Accordingly, each 
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Figure 34-- Undipped Projected Edges 




Figure 35- Edges Clipped in Y 



projected edge is assigned a fixed point X with indicator value to be set each time the 
edge crosses a scanline. 

10 X2<0 
n nn<X2<n, 0 < n < [2.57] 
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Utilization of for X-ordering reduces the infinite set {X!- «'<X<'H»} to the finite set 
{0,1,2, . , . , Rx}. 

If both Xj^> Rjj and X2> R^* the entire edge lies to the left of the viewing window 
and Ix^^X the edge, (a) If the right-side surface of the edge is a "flat-shaded" 
surface (i.e., the surface is planar and there are no point light sources in the world), the 
edge may be deleted from the set of projected atom edges since it would only be used 
to turn off a right-side surface of known constant brightness. That is, the left-side surface 
is not in the pyramid of vision and [1.15] is a constant, (b) If the right-side surface is 
"smooth-shaded," the edge must be kept to provide a brightness value and an X value 
for [1.15]. The edge's depth functions (for the surface on the left) are meanmgless. 

Similarly, if either Xj^> Rjj and X2< Rx or if Rx and 0<Xj^< Rx, the 
edge is used only for brightness data over the scanlines for which Ix"^X' 

If Xj> 0 and/or X2< 0, the edge must be modified. A new edge with X'-j^ =X'2 = 
0 = lx is created to represent the surface (which was started to the right of the viewing 
window by the original edge) along the right side of the viewmg window. The depth 
functions become 

D' = D. +Dx(0- X^) 
° ° ^ ^ [2.58] 

The original X and M are kept for the brightness algorithm [1.15]. Figures 36 and 
37 demonstrate the effects of horizontal clipping. Edges and lie to the left of the 
viewing window, the right-side surface of E-^ is "smooth -shaded" while the right-side 
surface of is "flat-shaded." Accordingly, E^^ is deleted from the final edge set while 
Ej is kept for possible brightness use. Edge E2 lies to the right of the viewing wirdow. 
A new edge with X = 0 and depth functions [2.58] is created to start the depth functions 
for the surface at the minimum X^0» the original X and M are still retained for possible 
brightness use. Edf-^ E^ extends both to the left and the right of the viewing window and 
generates two clipped vdgvb, Th< first generated Edge E'3 is at first useful for brightness 
only since Ix"^X ^^^^ ust ful for both brightness and surface depth. The second 
generated Edge E"3 provides an urate depth functions for its left-side surface at X=0, 
but retains the original edge slope for brightness use. 

Two clipping flags are also prc»vided in the ehpping section. Both flags are initially 
set to +1. The "left" flag, II, ii> cleared to 0 if a projected edge has X2>0. The "right" 



flag, Ip^, is cleared to 0 if a projected edge has X2<Rx- After all atom edges have been 
projected, the function I'=Il"**^R examined. Possible outcomes are: 

(1) 1=0: The projected atom surfaces are potentially visible. Add the 
projected atom edges to the y-ordered output data set for the visible surface processor. 

(2) 1=1: Either ^L"^ '^R'^^ projected atom lies either entirely 

to the left of the viewing window oi entirely to the right of the viewing window. Delete 
the entire projected atom edge set. 

(3) 1 = 2: Both Il = 1 and Ip^ = l indicating that all the atom edges failed to 
survive the y-clipping. The set of projected atom edges is empty. 

Window Clipping: Z 

Suppose an object edge in R^ has one end-point on the near (eye) side of the viewing 
window and the other end-point on the far side. There are singularities at z=0 in the pro- 
jection operators [2.21] and [2.22] , the projected edge slope [2.23] , and the depth func- 
tion [2.29] . Furthermore, if z^Zy^^, the range of the depth function is limited to [1,0). 
Therefore, projected objects shall be clipped to z^Z^^y much as they were clipped to an R^^ 
by Ry viewing area. 

The problem of clipping edges at z=Zy^ is somewhat similar to that of clipping 
edges on the right side of the viewing window at X = 0. The portion of the edge cor- 
responding to z> Z^^ may be projected into R and treated as a normal edge to be 
clipped in X and Y. The portion of the edge corresponding to z<Z^ may be ignored; 
however, this portion of the edge may generate a surface on its left side which extends into 
z> In order to start this surface and possibly to turn off a right-side surface, second 
and third edges are created. These new edges lie in the plane of the viewing window. 

Consider Figure 38, the intersection of plane ABCD and the plane of the viewing 
window, z-Z^y. The onginal plane is defined by the four edges: AB, EC, CD, and DA. 
It is readily apparent that: 

(1) Two different planes will either intersect along a straight-line 2 or will not 
intersect at all. Therefore, vsincr all atom edges defining a planar surface lie on that 
surface, all intersection>s of Ihoht* t'tJ^rs and the plane of the viewing window will he 
along some straight-lino 

(2) Since all atom edgps are te.-minated at points which are finite in R"^ and, 
since the edges form a dosed boundary fur the surface, there will either be no inter- 
sections or an f^ven number of inttirsections bulwi-en atom edges and the plane of the 
viewing window for vmh .surface, Sinci' the new <»dges formed between the points of 
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A^/ii//^ ^5 -Intersection of Plane ABCD and Viewing Window P?ane 

intersection between the edges and the vie^^ing window piano complete the boundarjr of 
the surface on the fa? side of the screen, all edges such as AB and BC which lie entirely 
on the ne^ side of the viewing window plane may be deleted, 

(3) If J is horizontal, all of these new edges lying in the viewing window 
plane are horizon ti\l aiid may be deleted; that is, projected edges of this type have no 
left or right surfaces. 

(4) If 2 is not horizontal, there is a unique y value for each point of intei- 
section and the new edges may be formed by starting at the greatest y value and 
connecting successive y intersections. Figure 39 demonstrates this procedure fo^ a 
planar surface which has a ''W" shape. Edges AB, CD, and EF are normal atom edges. 
Edge IJ lies entirely on the near side of the viewing window and is ignored. Edges GC, 
HE, lA, HD, IF, and GB intersect the viewing window plane z=Z.„„ at In lo I^ I.- 
and Ig. The points of intersection are y-ordered and new Edges Igig, I^I^, and I2I5 
complete the planar surface's boundaries. 
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Figure 39- Intersection of General Planar Surface With z = 2^^^ Plane 

The equation of line C can easily be determined from the surface normal, N, and a 
point, Pg, known to be on the surface. Let S denote the planar surface. Let P be any 
point on line 8, but !et P P2 as in Figure 40. 
Now since Pg^S and P€S, (P- is orthogonal to N and so 

(P-P2)"^"N = 0 

(x-X2)n^My-y2)ny+ {Z^^rh^""!^^ 

which yields 

' " (■^)'' ('^2-^y2^^(^2-2vw') [2.59] 



which has the form of a straight line, x = my + c. Since all points, of intersection between 
S and the plane of the viewing window satisfy [2.59] , Uu; slope of £ (m = -n^ln^) and 
the constant (0=^X2 + (Hy/n^jyg + (n.Jn^) (zg Z^^^i) need be calculated only once 
for each planar surface. Note also that m and c arc undefined only if 11^^ = 0. This occurs 
only if S and the plane of thp \iowiiig window interf-ecl along a horizontal line. Such 
horizontal cases are deleted s\s bt>foro. 
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Figure 40- Common Point to Both S and z = Z, 



Furthermore, each newly generated edge corresponding to surface S has depth 
functions d'^^l, d'g=0, and d'x=d^, where d^ is the incremental horizontal depth 
function for surface S. Figure 41 shows the data structure of intersecting surfaces and 
their corresponding y intercepts and edge parameters. Each time an edge that intersects 
the plane 'of the viewing window is discovered, a search is made for both the left- and 
right-side surfaces in the list of surfaces already discovered to intersect the plane of the 
viewing window. If a surface is found, the new y-intercept is inserted into the y-ordered 
intercept list for that particular surface. Also, the brightness at that position is calculated 
using [115] and retained with the y-intercept value. If a surface S is not found, then 
S, djj, m, and c are added to the surface list and its y-intercept and brightness at the 
point of intersection with the surface's first entries. 

After all the edges of an atom in have been projected, the surface data table 
is checked. If the table is empty, the projection processor goes on to the next 
([Tj] ,ATOMj) pair. If the table is not empty, the new edges are created for each 
intersecting surface passed thro gh the X-Y clipping unit. Only after this has been done 
does the projection processor proceed to the next ([Tj],ATOMj) pair. 
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Surface Data Table 



Intercept Data Table 



Surface 
Label 



m 



Pointer to 
Maximum 
Y Intercept. 



PNTR2 



0 



y4,b4 



PNTR 



0 



PNTR, 



Figure 41 - Y Intercept Data Structure: Yi > 72 > ¥3 > 74 



The convention for normal (Le., far-side) edges is that P2 produces a projected point 
above the projection of Pj^. This defined the left-right convention for the projected 
surfaces. The convention for edges that intersect the plane of the viewing window is that 
P2 is the end point on the far side of the viewing window. A second algorithm for 
determining the left-right nature of the surfaces for these edges is now needed. There are 
now three cases for projected edges: 

(1) Horizontal edges: Horizontal projected edges are still deleted, but they 
generate new surface edge intercepts at the projected y value. 

(2) Non-horizontal edges, dg<0: These edges 'O P'2 as the top projected 
point as well as P2 being the point of the far side of the v Jow. The lower projected 
point may be calculated by: 

Ay = d-dol'd^ 
Yl = 72 + Ay 
^1 = X2 + Ay " m 

where d^ is the depth function at r'2, the value of the dopth function at the plane of 
the viewing window is +1, and m is the slope of the projected edge» 



(3) Non-horizontal edges, dg> 0: These edges have P'^ as the actual upper 
point. Here 

Ay = (l-d^l/dj 
Vs = Vl 

Xj = 

Vl = yi + Ay 
x^ = x^ + Ay -m 

For this latter case, the left and right surface labels are switched as are (Bj^.Bg) and 
(63,3^), the brightness labels. 

The rule for determining whether a projected surface lies on the left or right side 
of a newly created surface edge depends solely on the d^ function for the surface. 
If d^<0, the surface depth gets greater as the surface is scanned from left to right 
(increasing x). Hence, the new edge is of the form 0:S. If d^>0, the surface depth 
decreases as it is scanned from left to right resulting in a new edge of the form S:0. 

Formulation of Atom Edge in r2 

A typical projected and clipped (in X, Y, and Z) edge and its relationship to the 
P-x ^y discrete resolution grid of the viewing window is shown in Figure 42. 

The depth function at the point (Xg.Yg) is D^. The incremental depth along the 
edge is given by Dg and the incremental depth horizontally to the left of the edge is 
given by D^. b^ and bg are the (log^) brightnesses of surface Sg on the left side of 
the edge and bg and b4 are the (log^^) brightnesses of surface on the right side of the 
edge, gj and g^ are the brightness gradients of Sg and along the edge such that 

b2 = bi +ge(Yi-Y2) 
^4 = b3+g^(Yi-Y2) 

The visible surface processor will determine surface depths and, hence, the visible surface 
only at discrete positions on the resolution raster. This limits the number of times that 
surface depths need be calculated and compared to a maximum of • R„ times if the 

y 

depth is to be calculated at every discrete position. 

Consider surface Sg as defined by the edge of Figure 42. The first scanline for 
which the edge becomes valid starting surface Sg to its left is Scanline Y-j. Subsequent 
scanlines for which the edge starts surface Sg ai-e;!i; 5=^+1, j+2, . . . , k- l,k. Since the 
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Figure 42- Projected Edge on Discrete Raster 

edge starts Sg only over this more limited Y range, the edges are equivalent to a shorter 
edge that starts exactly on Scanline Y=j and ends exactly on Scanline Y=K. This 
shorter edge has adjusted parameters given by 



\ = X2 + M(Yj-Y2) 

D'o = Do ^Ds<Yt-Y2) [2.60] 
b', = b, +gg(Yt-Y2) 
b'3 = b3 + g^(Yj-Y2) 

The resulting edge is shown in Figure 43. With the shortened edge of Figure 37, the depth 
of the surface along the edge at any scanline n scanlines from the top scanline for which 
the edge is valid is given by D=D'o + nDg. In particular, a /ule for updating the depth 



function from one scanline to the next is now given by 



(2.61] 
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Furthermore, the position of the edge at any scanlinu n scanlines from the top scanline for 
which the edge is valid is given by X = X^ + nM, In particular, a rule for updating the 
horizontal edge pdsition from one scanline to the next is now given by 

- + M [2.62] 

The first scanline for which the edge is valid is Scanline j where j is the first integer 
greater than or equal to That is, j l<Y2<j. The last scanline for which the edge 
is valid is Scanline k where k is the first integer strictly less than Yj^. That is, k<Yj^<k+l. 
This definition assures that a surface will never appear twice on the same line due to two 
edges terminal ng at a vertex point that, by chance, is projected exactly onto a scanline. 
This is demonstrated in Figure 14. This procedure of clipping the tops and bottoms of 
projected edges exactly at scanlines is eabil> accumplibhed by defining a height function 
H for the edge 

H- (ib- ib^l<Yi<ih [2.63] 

7/ it-i<Y2"it 
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Figure 44 -Pro'iected Edges Ending on Discrete Scanlines 
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If H^O, edges are clipped as in Figure 44. If H^O, the entire projected edge lies between 
successive scanlines. Since these edges fail to intersect any scanlines, they are deleted from 
the projected atom edge set. When an edge first becomes valid at Scanline Jt, its height H 
is as given by [2.63] . Then, for each successive scanline, the remaining height is given 
by H-1. So the rule for updating the height from one scanline to the next is given by 

H = H-1 [2.641 

H=0 corresponds to Scanline Jb and indicates the last valid scanline for the edge. 

The last parameters to be updated from one scanline to the next are the brightness 
values. These are updated by 

[2.65] 



[2.66] 



b3 = b3 + 

over the range vi^here H^O. 

Thus, once the top of a projected edge is clipped to a Scanline Yt=j|^ via [2.60] 
and a height function H is defined via [2.63] , the edge is fully defined by the set 
of parameters 

where I^ is the valid X indicator function defined by [2.57] for each scanline that the 
edge crosses. Also to be updated at each valid scanline are H, D^, bj, and b^. These 
updates require only addition operations between successivi' ocanlines. 
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One further modification has been arbitrarily made. Atom surfaces for worlds con- 
taining no point light sources are typically flat-shaded. These surfaces have a constant 
brightness within their edge perimeter. Accordingly, edges with such surfaces on the right 
side h? e b3=b4 and gr=0. These edges are thus used only to end the right-side surface of 
know., brightness. Therefore, the r2 edge convention is taken to be 
(Yt,H,Ix,Xt,M,S£,-Sj,Do,Ds,Dx,bi,gg,c,Is) where Ig is a flag indicating whether or not 
surface Sg is smooth -shaded or flat-shaded. For cases where is a smooth-shaded 
surface requiring a brightness function b2,gf different from bi,gg, a second edge is 
formed with parameter (Yt,H,Ix,Xt,M,-SB,Sr,Do,Ds,Dx,b3,gr,c,Is). 

Y-Ordered Output Edge Set 

Projected atom edges that survive all the clipping processing are added to an output 
edge set. This output edge set is ordered according to the Yj-'s of the edges. The two 
tasks involved in adding a projected atom's edges to this set are (a) maintaining the 
Y-order of the Y^'s and (b) relabeling the surface labels in the projected edges to 
maintain uniqueness. 

The Y-order is maintained through the use of a one-pass bucket sort technique. 

Two tables are employed. Table Y has Ry positions-one for each of the R„ scanlines. 

•' y 

Table E has E^ax locations, where E^ax >s the maximum number of projected edges 
allowed. Both tables are initially cleared. As each edge is added to the output set, Table 
Y's Yjth position (Y^ is the starting scanline for the edge to be added) is checked for 
data. If it is null, the edge being added is the first occurrence of an edge with that 
starting scanline value, and the edge's position q in the output table is placed at 
position Yt in Table Y. If position Yj is not empty but contains some number p, then 
the pth edge in the output buffer also has Yj as its starting scanline. Accordingly, the 
null in the qth position of Table E is replaced with a p and the p in the Y^^^ 
position in Table Y is replaced with a q. After all edges from all the projected atoms 
have been entered into the output table. Table Y is scanned for data. If the j*^ position 
of Table Y contains a k, the k*^ edge in the output table is r-ad out. Then, position k 
in Table E is checked for additional edges starting on scanline j. If position k in Table 
E is null, the reading of Table Y is resumed at position j+1. If position k in Table E is 
not null but contains C, Edge 2 in the output table is read out and position £ in Table E 
is checked as before for additional edges starting on scanline j. 

Unique surface labeling is achieved as the edges are added to the output edge table. 
Initially an offset number K is set to one. Then, as each atom's projected edges are 
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added to the output table, +K is added to positive surface labels and -X is added to 
negative surface labels. 0 labels remain 0. After each atom's edges have been added to the 
itput table, K is incremented by Ns, where Ns is the highest surface label encountered 
in the atom last processed. For example, suppose the world consists of two cubes, a 
cube having six surfaces which are labeled in the atom library 1 through 6. Then the 
edges for the first cube will enter the output table with surface labels in the range 
2 through 7 while the edges for the second cube will enter with surface labels 8 through 13. 

Two special edges are always added to the output table for each display. Each has 
H=Ry, M=0, Do=l, D,.=D, Ds=0. The first of these edges has Xt=0 and surface 
labeling 0:1; this edge serves to initialize the display at the rig.it side of the viewing 
window. The second edge has X^=R^ and surface labeling 1:0; this edge serves to mark 
the end of the display at the left side of the viewing window. It is the presence of these 
edges in every display that initiates the value of the surface label offset number K at one 
rather than at zero. 
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SecJion 3 
VISIBLE SURFACE PROCESSOR 



INT RODUCTIOIM 

This section provides a detailed description of the visible surface processor of the 
CHARGE image generator. The visible surface processor performs two functions: 

(1) It solves the visible surface problem at each of the by Ry resolution 
points on the observer's viewing window. This process generates the video 
image of the world as viewed by the observer. 

(2) It encodes the video image into a data compressed format, allowing it to 
be buffered locally at individual CHARGE terminals at a low cost. 

Since all displays are locally buffered, the image generator and the communication link 
to the terminals may be time shared. 



INPUT DATA 

Input data to the visible surface processor is the Y-ord3red output edge set of the 
projection processor representing potentially visible projected atom surfaces* Each of 
these projected edges is a 14-tuple of the form 

Ep = (Yt,HJj^,Xt,M,S8,S,,Do,D5,Dx,b,g,cJs) [3.1] 

where 



yt 


' First scanline for which Ep is valid 


H 


: Height of projected edge Ep 


Ix 


Valid X indicator function defined by [2.57] 


Xt 


Horizontal position at Scanline Yj. 


M 


Slope of Ep (M = AX/AY) 


Sg 


Label of left-side surface 


Sr 


Label of right-side surface 


Do 


Value of depth function at (Xt,Y|;) 


Ds : 


Incremental change to depth function along edge 


Dx : 


Incremental change to depth function horizontally 


b • 


log brightness at (Xt.Yt) 
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g : log brightness gradient along edge 

c : log color of Sg 

Ig : Indicates if Sg is smooth-shaded 
A detailed explanation of these parameters is given in Section 2, Formulation of Atom 
Edge in r2 (pg. 72), 

OUTPUT DATA 

Output data to CHARGE terminals consist of Y -ordered (top tc bottom) blocks of 
X-ordered (left to nght) visible edges. Each of Ihrse visible edges is an 8-tuple of the form 

Ev - (Yj,H.c.X^.Mb,gJ^) [3,2] 

where 

: First scanline for which liy is valid 
H : Height of visible edge Ey 
c : log of color to right of visible edge Ey 

• Horizontal position of Ey on Scanline Y^ 
M : Slope of Ey (i\I=AX/aY) 
b : log of brightness at (Xt»Y|.) 
g : log of brightness gradient along Ey 
Ig : Indicates if right side is to be smooth-shaded 
Each such visible edge controls the color and the brightness of an area to its right; 
hence, the X-ordered set of edges within a block controls that area of the TV screen 
for which block is vidid. The first edge within a block i;> a specially encoded **marker 
edge" which ser\'es to both identify the start of the block and to sot the Y-limits over 
which the block is valid. Furthermore, special ''null edges*' may be located anywhere 
within the block. The number of edges within a block is bounded by 

Nb - N (256- iAY/401) [3.3] 

where 

N : Numl^er of parallel memory hanks assigned to the CHARGE terminal. 

AY : Y-range of the block 

[AY/40] : Truncated integer 
A detailed description of the codmt^ of the eight vislblr cd^o parameters into an Tfi-bit 
word and the CHARGE terminal memory ^tru( turc is found in a separate report.* 
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■ OVERVIEW OF SCANLINE OPERATION 

The operation of the visible .surfacf proct-.v-oi- i.-, on a sun esbive scanline basis whereby 
each scanline represcrt:. a di^cn-te Y rtssolutiua \ ilue of tb" i.twrver's viemng window 
(Figure 32). This operation is as follows: 

(1) Starting at the top of tho observer's viewing window (Y=l), the visible 
surface processor partitions each t.,. online i.i the X interval [Rx.i; nto a set of mutually 
exclusive, collectively exhaustive, subinterval.-> of the form (X,+j,X,] such that either no 
surface or exactly one surface is \x<M,' m e.uh subinterval. This is accoinphshed by deter- 
mining the nearest ■.>uri;ice for ;J1 di^dvte wdu^s of X in the interval [R^,!]. 

(2) As it partitions successive stanlines, the visible surface processor creates 
and maintains an X-ordeved set of straight bne visible edges that affect the same scanline 
partitioning as in (li. As the scanline pariit,oninf: chaiv^es from one scanline to the next, 
some edges in this set remain v.ilid and may t-Monded, other edges cease to be valid 
and are deleted from the < urrcnt .,et. and still othi>r vi^ibk' edges must be created and 
added to the current visibU- edge set to realize the . urrent scanline partitioning. Visible 
edges deleted from the current xiAhli rd?.- set remain as X-urdered members '^f the output 
visible edge block whose size i.^ I)ounded hy [DM]. 

(3) Step. (1) and (2) ;ire repealed uniil Y-Ry. At this point, the final block of 
visible edges is coinplete and ; .Klitional bloc ks fill. d with marker edges and null edges are 
created to fill up the CHAUCJC terminars memory- lianks and to provide null data for 
those scanaiies c-Tespondn- 1,. vertua! retrace iime of the CHARGE terminal's monitor. 

Figure 45 demonstrates for a stnul^ Scanhnc- Y - j the data flow from the input data 
set of Y-ordered projected ed.;.«s to the r atput blo.di of X-ordered visible edges. At the 
start of Scanline Y-i, any cdv'. s ir, tl^.. mpm data set that !)ccome valid at Y=i are added 
to the set of projected etige. that re adJ.>d on preceding scanhn(>s but which remain 
voJid (due to their height v.-.hu-s, 11, . Con.panif.n to this current projected edge buffer is 
a pointer list that roaint uns a uuhf t-, l. fi X-ord-nnf,' <.f tho projc.-iod -dse set at the 
resolution of Rj. dis. rete X wI<'T\-A-. per .,:m\mf\ Ki. h edne m the current projected 
edge set points to a next rd-'^ v.ho^e X valu*' hes iMthcr v.-lbin the same X resolution 
interval or in a higher X resdiition n..'.!. 

Thf visible surf-j. e aI"onfhn ^.r x th.-,- .,!.<(,. m X ta-df-r from njiht to lej't, 

genenlmg hon/ontai i;it<-p;;;-. r.iiud (i,h.'t (i, uiri i,r ih.- ,.,ami> -.■iirface is visible. 

Each edge as it is encountn ed tnes to .-nd output atleinf'is f..r (he surface on its right 
side and to .start output attfrnpf - for \hv ^.irf.... on if. !.-f( side T!..- \isihle surface 
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Figure 43 - Visible Surface Processor Functiono! Overview 
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algorithm utiUzes the surface depth parameters included in the projected edge definition 
to determine that surface whose distance from the observers eye is minimum. When the 
niimmum depth surface changes, the visible surface algorithm forms a p.^xtition between 
Uie two minimum surface intervals; the attributes that it assigns to the partitioning 
boundary are those that are necessary to uniquely match a visible edge partition. 

At the end of the scanline (X=R^) the partitioning is complete and the set of current 
projected edges may be updated for the next scanline (Y=i4.i). This updating process for 
each edge involves: 

(1) H=H-1. If H<0, the edge is not valid for Y=i+1 and may be deleted from 
the current valid projected edge set. 

(2) Do=Do+Ds. This updates the'depth function at the top of th^ e^'^". 

(3) X=X-!-M. This updates the X position of the top of the edge. 

(4) Ix, the X indicator function as defined by [2.57] is updated by the new 
X value. 

(5) The X-ordered list is updated by the value, if necessary. This crossover 
of edges is a relatively lov/ frequency event and does not require a complete 
X-ordering of the entire current projected edge set. 

In essence, the update procedure "shortens" each edge by clipping off the top of an edge 
and redefining the new top in X,Y, and depth. 

Analogous to the set of current projected edges at Y=i and its companion X-order 
list are the set of current visible edges at Y=i-l. When the actual video scanline partitioning 
has been completed, an attempt is made to match each successive actual video partition with 
a visible edge extended to scanline i.The order of this matching attempt is governed by the 
X-ordering of the preceding scanlinc's visible edges. As visible edges are matched, they 
are deleted from the Y=j-1 ordering list. 

If a matched visible edge is not the first edge in the remaining Y=i- 1 order list, those 
edges preceding it in the list are deemed nof to be extendable onto scanline Y=i (since 
visible edges may never cross) and axe ended at scanline Y=i -1. If an actual scanline 
partition has no corresponding visible edge, it is deemed to be "new" and a visible edge 
is started for it. This new or matched visible edge is next compjired in X with the pre- 
ceding edge to determine if it is "too dose." Visilile edges surviving this comparison 
generate a visible output table for scanhne Y=,. At the end of the scanline, this table is 
used to generate an X-order list for the current visible edge data set and to insert into 
proper X-order any new visible edges added to the visible edgo block. 
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The visible edge block is ended and sent (read out according to its X-order list) 
whenever the number edges in it plub the number of actual video partitions exceeds 
N3 as defined by [3.3]. This guarantees that each scanline is uniquely defined by only 
one visible edge block. As a visible edge block is read out, each of its positions except 
for those corresponding to udgub that still appear in the current visible edge set is freed 
up. Hence, the initial X-ordering for the next visible edge block is that of those edges 
extending down from the preceding edge block. Thus, if the image is as in Figure 46 
and blocking occurs at Y=Yj^, the corresponding visible edge blocks are those of 
Figure 47 and 48. Note that while edges do not extend down into subsequent blocks, any 
edge may have its origin defined in the Y-range of any preceding block. 

VISIBLE SURFACE ALGORITHM 

Scanline operation of the visible surface algorithm is shown in Figure 49. Initially, 
the output table, the pending delayed remove table, the change table, and the pending 
brightness table are cleared. X, the horizontal position register, is set to zero and the 
first "beam box" is loaded with surface and the others are cleared. This first beam 
box will always be cleared by the 0.1 edge at X=l that is present in every image to turn 
on the appropriate background color and brightness (see Section 2, y-Ordered Output 
Edge Set, pg. 77). Also initialized at the start of each scanline are 12 memory registers 
that are used to determine for each X whether or not a now video partition should be 
added to the output table: 



(1) 


ICNAM 


: Surface with minimum depth at X-i 


(2) 


ILNAM 


: Surface with minimum depth at X=i 


(3) 


ICPNT 


: Beam box containing ICNAM 


(4) 


ILPNT 


: Beam box containing ILNAM 


(5) 


DC 


: Depth of ICNAM at X=i 


(6) 


DL 


: Depth of ILNAM at X=^i-1 


(7) 


ICSM 


: Smooth-shading bit for ICNAM 


(8) 


ILSM 


: Smooth-shading bit for ILNAM 


(9) 


ICC 


: Color oflCNAM 


(10) 


ILC 


: Color of ILNAM 


(11) 


BC 


: Brightness of ICNAM 


(12) 


BL 


: Brightnc-^s of ILNAM 



At the start of each soanlino ICNAM=1. ICPN1--1, DC==1, anrl ICSM-0. 
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Edqes Defining Image 
Figure 46 ~ Edges Defining Image 




Range of 
Y to be 
Used for 
Block No. 1 



Edges m Block 1 
Figure 47 - Edges in Block 1 




Edg«s JO Block 2 
Figure 48 Edges in Block 2 
87 



Horizontal 
Skip Ahead 




NO 



® 


Add Edge Data 


to Change Table 


® 


Add and Remove 


Beam Boxes 




According to 




Changft Table 




Update Delayed 




Remove Table 



Update Beam 
Box Depths 
Getting New 
{CurremJ 




Add Video 
Partition 
to Output 
Table 



Fi:)'ire 49 ~ Visible Surface Algorithm 



After initialization, a processing loop transforms the input set of X-ordtred pro- 
jected edges into an X-ordered output set of visible surface partitions (the visible surface 
between partitions is t\at surface whose depth is minimum o\er Lhat X interval between 
partitions). This processing loop operates only at the discrete values X=ie(i,2, . . R^}, 
and the operation is in one of two modes depending upon whether or not there are 
projected edges in the interval [i.i- 1). If edges are present in that interval, those beam 
boxes containing the edges' right-side surfaces are unloaded and reloaded with the label 
and depth information of the edges' left-side surface, all of the beam boxes are then 
updated to X=i+1 and the surface with minimum depth at X=i+1 is determined. If the 
next edges are present at X=j (j>i), all of the surfaces that are potentially visible at 
X=i+1 are potentially visible at X=j 1, hence, a horizontal skip-ahead attempt may be 
made by updating each beam box by (j-i-1) units. If the same surface visible at X=i is 
visible at X=j-1, the skip-ahead is known to be valid since all surfaces are planar. 

The format of a beam box is shown in Figure 50. The depth of the surface in the 
perspective domain is given by parameter D(1>D>0). is the incremental change 
in the perspective depth in the horizontal direction and DSAVE is a register used to save 
the value of D during a horizontal skip-ahead atcempt, D is restored from DSAVE should 
the skip-ahead attempt prove unsuccessful. LABEL identifies the surface, ISM is a bit that 
is turned on if the surface is smooth-shaded, and PNTR points to the edge in the 
X-ordered current projected edge set that was used to turn on the beam box. The number 
of beam boxes loaded with surfaces at X^i is tho number of overlapping surfaces at X=i. 

LABEL 
ISM 
PNTR 
D 

DSAVE 

Figure 50 - Beam 
Box Format 

Each beam box is loaded in such a way that its depth parameter D represents the 
surface perspective depth precisely at the discrete values of X along the scanline. This 
assures that all surface depth parameters will bt consistently compared. Consider Figure 51. 
Edge E crosses the scanhne J at i>X^>i 1, causing the edge to be encountered in the i^h 
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Figure 51 - Edge E in ith X Range 



X range. The depth function D is to be loaded such that it will be a correct measure of 
the left surface depth at X=i,i+1. . . . after updating. This is easily accomplished by setting 

D = Do+lrX^lD^-Dj^ 

= Do+(i-Xt"l)Dx [3.4] 
That is, the beam box is loaded as if the left-side surface existed at X=i- 1; then, as the 
beam boxes are incremented forward 1 unit by adding to D, all the b^xes refer to 
depths at X=i. 

Often more than one edge may reference the san^e surface within an X resolution 
interval. This situation occurs when edges converge at a vertex point of an atom-- 
particularly if the projected atom size is small and of the same order of magnitude as 
the rastering (W^/Rx). To handle cases such as the.se, all of the edges within the same X 
interval are scanned before any changes are made to the beam boxes. This scanning pro- 
cedure makes entries in the Change Table, each entry in the Change Table has the form 
shown in Figure 52. 



LABEL 
Net Change 
PNTR1 
PNTR2 



Figure 52 - Change 
Table Entry 
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As each surface is encountered in edge definitions, a running count is kept of the 
number of times that it appeared on the lefl side of an edge (add) minus the number of 
times that it appeared on the right side of an edge (remove). Then 

Net change = # adds - # removes [3 5] 

If net change is positive, a beam box for tae surface is added (if one does not already 
exist). If net change is negative, the corresponding surface beam box is removed (if it 
exists). If net change is zero, no attempt is made to either add or remove a beam box for 
that surface. If an attempt is made to remove a beam box that does not exist, that attempt 
will be remembered in a "delayed remove" table and will be used to negate any single 
attempt at a successive X location to add a beam box for that surface. These delayed 
removes are caused by slightly inaccurate (-usually very large) projected edge slopes as in 
Figure 53. PNTRl points to the edge (at X=i) in the current projected edge set with 
LABEL on the left side and (a) ILNAM or (b) minimum surface label on the right side. 
PNTR2 points to the edge (at X=i) in the current projected edge set with LABEL on the 
right side and a minimum surface label on the left side. Usually either PNTRl is zero 
(remove only) or PNTR2 is zero (add only). 



X i+1 i i-i 



M 
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Figure 53 - Delayed Removal of 
Surface at X = i + 1 on Scanline j 
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Once all of the changes have been made to the beam boxes and the beam boxes have 
been updated by one X resolution level, a comparison of the depth parameters is made 
to determine which beam box contmns the surface with minimum depth. If two surfaces 
have the same minimum value of D, selection of the surface with the minimum surface 
label guarantees consistency in resolving such conflicts. Selection of a beam box sets 
ICNAM, ICPNT, and DC. The pointer contained in the beam box is then used to obtain 
values for ICSM, ICC, and BC. 

At this point, the test for change in output is made; ICNAM and ILNAM are com- 
pared; if ICNAM=ILNAM, the surface visible at X=i-1 remains visible and no changes to 
the output set are required. If, however, ICNAM^^ ILNAM, a change in visible surface has 
occurred in the ith X interval and a new partitition must be added to the output video 
partition table. The partition between ICNAM and ILNAM is a straight line edge in the 
plane of the vie\ving window and may be caused by one of three situations: 

(i) ICNAM may have been added by an edge of the form ICNAM :S in the 
i^^* X range, where S is any surface label. Typically, however, either S=ILNAM corres- 
ponding to an edge between surfaces on the nearest atom, or S is negative corresponding 
to the rightmost edge of a new atom that lies closer to the eye than that atom con- 
taining surface ILNAM. 

(ii) ILNAM may have been removed by an edge of the form S:ILNAM in the 
i^^ X range, where S is any surface label. Typically, either S=ICNAM corresponding to an 
edge between surfaces on the nearest atom, or S is negative corresponding to the left- 
most edge of the last visible atom that now ceases to be visible. 



In Figure 54, scanline j is partitioned into seven regions by three surfaces. S2 lies 
in front of S]^, and S3 intersects Si in a line passing through point c. Each partition is 
also labeled as to being cau.sed by Cases (i), (ii), or (iii). 

The partitioning is: 



(iii) There may be an intersection between ICNAM and ILNAM. 



(1) 
(2) 
(3) 
(4) 
(5) 
(6) 
(7) 



(a,l] 
(b,a] 
(c,b] 
(d,c] 
(e,d] 
(f,e] 

tRx,n 



ICNAM=Si 
ICNAM=S3 
ICNAM=Si 
ICNaM=S2 
ICNAM=Si 
ICNAM=0 



ICNAM^O 
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where 0 denotes "background only*', that is, no projected atom is visible and the video is 
the background color and brightness. The partitions at X=l and X^R^ are Cases (ii) and 
(i), and are due to the 0:1 and 1:0 edges present in all displays. 

The format for a video partition L shown in Figure 55. This format contains all of 
the information required to match an extended visible edge or to start a visible edge 
should the partition fail to match an existing visible edge. Sj^ and are left and right 
surface labels for the partition. Si^^^ is set to ILNAM and COLOR=ILC, thereby referring 
to the surface to the right of the partition. TYPE=2 ICSM+ILSM, thereby indicating which 
of the following smooth shading transitions has occurred: 
TYPE = 0 = 00 -> flat : flat 

^ 1 = 01 -> flat : smooth 

2 = 10 smooth : flat 

3 = 11 -> smooth : smooth 

If the partition is due to a projected atom edge as in Case (i) or (ii), one or both of 
the NEW/LAST bits may be set. The NEW bit is set if the edge in the current projected 
edge set has become valid for the first time; that is, the edge has just joined the current 
projected edge set. The NEW bit set causes the partition to generate a new visible edge 
since any match with an existing visible edge would obviously be erroneous. The LAST 
bit is set if the edge in the current projected edge set has 0, that is, the edge will be 
deleted from the current projected edge set at the conclusion of the present scanline» The 
LAST bit set causes its corresponding visible edge (matched or new) to be terminated on 
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^last 
COLOR 
TYPE 



NEW/LAST 



Figure 55 —Video 
Partition Format 



the present scanline since any extension downward to succeeding scanlines would oDviously 
be erroneous. Xp.Mp, and Bp reflect the position, slope, and brightness (to the right) of 
the video partition. 

All of the data required for a video partition are available at the instant it is created 
except, perhaps, for the value of bp. If the partition is due to a projected atom edge, 
then S£,Sp, NEW/LAST, and Mp are obtained from that edge while Sjagt, COLOR, TYPE, 
and Xp are obtained from ILNAM, ILC, ICSM, ILSM, and X. If the partition is due to an 
intersection, then Sg=ICNAM, Sr=ILNAM, NEVV/LAST=00, and Mp is calculated from the 
edge data of the edges pointed at by beam box #ILPNT and beam box #ICPN1 . Consider 
Figure 56. Edge El with Slope M2 loads a beam box with SI and Edge E2 with Slope M2 
loads a beam box with S2. At X=i, ICNAM=S1 and ILNAM=S2. At X=Xp on Scanline j, both 

51 and 82 have a perspective depth of D=D'. The dashed line of intersection between 81 and 

52 represents the locus of points on 81 and 82 having the same perspective depth. 

In particular, at Y=j+1 

D(Si) = D(S2) 
D'+Dy^+Mp.D^^ = D'+Dy2+Mp.Dx2 
Dy,+Mp.Dx^ = Dy2+Mp.Dx2 
But, since 8^ and 82 are planar, Dg. = Dy.+Mj.Dxj and so 

or 

«Ds2-'^2-Dx2)-(Ds,-Mi.Dx,) 



Mp 



[3.6] 
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Figure 56 - Intersection of Visible Surfaces 



If surface ILNAM is not a smooth-shaded surface, bp=BL since the brightness at any 
point on a flat-shaded surface is that of the surface when the beam box was loaded by a 
right-side edge. If ILNAM is a smooth-shaded surface, and there exists an edge of the 
form S:ILNAM within the interval (i,i-l] where i^^Rx, then the brightness at X=Xd is 
given by the brightness of the edge; if, however, i=Rx, then the edges' horizontal position 
discretely limited to by the indicator function I^ (given by 2.57) for X-ordering 
reasons, may not be R^ in which case bp is given by 

bp =br+iV^ (bjrbr) = br+ {Xp-Xy)g [3.7] 
(Xu-Xr) 

where the subscripts 8 and r designate the left and right edge parameters and Xj and X^ 
may lie anywhere on the scanline subject only to Xg>Xj. 

If no edge of the form S:ILNAM exists within the i^h x interval, the value of b^ 
cannot yet be ascertained. In these cases, the brightness is deemed to be pending and an 
entry is made in the "pending brightness table," which has the form shown in Figure 57 
for surface Si with three pending brightnesses at positions Oi, O2, and O3 in the output 
video partition table. When an edge of the form S:Sj is encountered yielding a remove 
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beam box Sj, the Sj of the edge is matched against the Sj in the pending brightness table 
indicating that one or more pending brightnesses must be calculated. The horizontal 
brightness gradient is first calculated using the bnghtness and X position data of the two 
edges that load and unload the beam box 

^g"br [3.8] 



9 = 



Once this gradient has been calculated, each of the bp's in the pending brightness table 
may be calculated by [3.7]. Referring back to surface SI in Figure 54, the brightnesses 
to the right for the partitions at X=b and X=d are calculated at X-f if surface is a 
smooth-shaded surface. 



Sl 


Pi 


S2 


P2 






Si 


Pi 






Sn 


Pn 


0 









02 



O3 



Figure 57 - Pending.Brightness Table Form 
k " 

If ICNAM is a smooth-shaded left-side surface, and there is no edge of the form 
ICNAM:ILNAM in the i^^ x interval, the image generator must form two partitions. The 
first partition is at X=Xp and affects the desired right-side brightness and color. Further- 
more, TYPE=ILSM controls smooth shading across the right-sidi' surface. The second 
partition is formed at X-Xp+2 (minimum separation for the CHARGE termmal); this 
second edge has TYPE=2, COLOR^ICC, and is used to provide a brightness for the 
CHARGE terminal to smooth-shade toward as it decodes the left-side surface. 
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) 

The addition of this edge at X=Xp+2 is always tentative and will be deleted should 
a partition occur within [Xp+2,Xp) with a right-side surface which is not equal to 
the left-side surface Sg of the double partition at Xp. If Sr=S£, the ne^^/^artition is 
deleted and the old partition remains valid. Figures 58 and 59 demonstrate hc^i of these 
cases. Surfaces Sq, Sj, and S2 are surfaces on a cylinder atom; Sq forms one of the flat 
cylinder ends whUe and S2 are t vo smoothly shaded facets of the cylinders round outer 
surface. Surface S3 lies on a olane nearer to the eye of the observer than the cylinder. 

On Scanline i, a double edge is created to account for the smooth-to-flat S^iSq visible 
surface switch. Farther down the scanline there is a S2:Si partition (no brightness dis- 
continuity), and still farther down the scanline is a S3:S2 partition. On scanline i+1, 
the double partition is initially formed to affect the Si:So visible surface switch; however, 
before X can increment by 2, the S2:Si visible surface switch occurs. Since surface Sj 
appears in both places, this second partition may be safely deleted since the brightness 
function defined over the curved cylinder is continuous. 

At Scanline k, however, the S3:-S4 (S4 is the back surface of the plane) partition 
occurs within two X resolution intervals of the S2:So partition. In this case, to allow the 
TYPE=2 S2:So partition to survive would be to erroneously move it across a brightness/ 
color discontinuity. Figure 59 demonstrates the final correct video partitioning of this 
region of interest with arrows indicating smooth shading. Note that on scanline k the 
S3:-S4 edge incorrectly shades toward the brightness of Sq, the flat end of the cylmde., 
rather than toward Sj, the rounded surface. 
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Figure 58 - Projected Plane and Cylinder 
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Ftgurf^ 59 - Close-up View of Plane and Cylinder 

This error will be unnoticed by the eye since either: 

(i) The distance between the two partitions is less than two X resolution 
intervals (out of 1600). 

(ii) The S3:-S4 partition may be shortened to scanline k-1 by the matching 
algorithm, thereby causing flat surface S3 to "jut** over by this small distance to Sq on 
scanline k. 

If a video surface switch is due to an edge, as in Case (i) or (ii), care must be taken 
to select the correct edge. Consider, for example, Figure 60, Surfaces Si,S2, and S3 all lie 
on the same atom. At X=i, there is a visible surface switching from Si to S3, Within the 
i^'^ X interval are two edges: of the form S2:Si, and E2 of the form S3:S2. Edge Ei 
is Case (ii). Edge E2 is Case (i). The problem is to select the parameters of one of the 
edges for the partition, 

A second ambiguous case is that of Figure 61 in which several edges reference Sj^ as 
the left-side surface (and S2 as a rignt-side surface). 

A third case is that of Figure 62. In the first case, the new visible surface, S^, actually 
lies in front of S2 while in the second case, S2 is in front of S^. 

These ambiguous situations may be consistently resolved by the processing steps of 
Figure 63 which utilize ICNAM, ILNAM, DC, and DL, NBEAM contains the number of 
beam boxes vying for output. 
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Figure 62 ~ Choice of Edges at X = i 
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Figure 55 - Setting Output Mode 
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The intersection mode is a default mode that is entered if either there are no edges 
in the i'h x resolution interval or else none of the edges in the 1*^ resolution interval 
reference ICNAM or ILNAM. The partition has Sj=ICNAM and Sr=ILNAM, and Slope M 
as calculated by [3.6]. One partition is created if ICNAM if flat-shaded; if ICNAM is 
smooth-shaded, a second partition is created at X=i+2. If either surface is smooth-shaded, 
that surface is added to the pending brightness table. 

The add and the remove modes are quite similar. If both ICNAM and ILNAM are 
referenced by edges, DC and DL are compared to determine which of the two cases of 
Figure 62 has occurred. Regardless of whether the add or the remove mode has been 
entered, a ICNAM:ILKAM edge always has priority: If such an edge is found, IMATCH is 
set to point to it and its parameters are used in forming the partition. Otherwise, if the 
add mode is entered, the edge whose parameters are used is that edge with the smallest 
right-side surface label. If the remove mode is entered, the edge whose parameters are 
used is that edge with the smallest left-side surface label. Any error made on this scanline 
will be corrected on subsequent scanlines as the edges diverge into different discrete X 
intervals. If the edges do not diverge, the error is not significant enoujih to warrant 
correction. 

Since the number of projected edges per scanlines is typically much less than Rx, 
the number of X resolution intervals, and since the occurrence of intersections involving 
front surfaces is a low-frequency event, the visible surface algorithm of the image g( lerator 
has been given the ability to horizontally skip-ahead over regions void of projected edges. 
The initial skip-ahead attempt is to the X interval preceding the next edge. That is, if 
the next edge occurs at X=i+AX, then 



where K is the number of surfaces vying for output. If ICNPNT=ILPNT, the skip-ahead 
was successful and the visible surface algorithm may now go back to a one unit mode 
after loading and unloading the beam boxes in accordance with the edges at X=i+AX. 
If ICPNT^ILPNT, there is at least one intersection in the interval [i+AX,i] and the 
maximum likelihood estimate is that there is exactly one intersection. The point of mter- 
section between the surface visible at X=i and the surface visible at X=i+AX may be 
calculated by equating their depths: 



Dk = Dk+(AX-1)D^. k=1 K 



[3.9] 



DSAVE^+D^^.DELTAX = DSAVE2+Dx2-DELTAX 



DELTAX = DSAVEg-DSAVEi 



[3.10] 
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where DSAVE is the register in each beam box remembering each depth D before 
skipping ahead. The next skip ahead increment is [DELTAX], the first integer less than 
DELTAX, and the process is repeated. Choosing [DELTAX] in this manner guarantees 
that the intersection will be caught between both ends of a one interval X increment. 
Figure 64 demonstrates how the intersection between Surface Sj, which is visible at 
X=i, and Surface S2, which is first visible at X=[Xk'] + i, is detected in three steps. 

The first skip-ahead attempt is made at X=j- 1, the first in.ervai before projected 
Edge E4. The visible surface is S3 necessitating the calculation of X^, the point of inter- 
section between surfaces Sj and S3, The boxes are next updated for X=[Xi^] ; here, the 
visible surface is still not Sj but S2. The point Xj^' is then calculated and the beam 
boxes are updated for X=[Xk'] ; here, Sj is still the visible surface and an additional 
one-unit increment obtains the surface ordering on the other side of the S2:Si 
intersection. 




Eye 

Figure 64 - Three Step Skip-Ahead 
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As each point of intersection Xj^ is calculated, it is checked against the 
maximum point X = j - 1 to assure that numerical inexactness (particularly with nearly 
coplanar surfaces) has not led to an Xj^ value that lies outside the permissible scan- 
line range. 

VISIBLE EDGE ENCODER 

Once the X-ordered set of video partitions for scanline Y = j has been form«3d, it is 
sent to the visible edge rncoder to be compared against an X-ordered set of visible edges 
that effect similar partitions on previous scanlines. Each visible edge that is matched by a 
video partition and is not "too close" to the first visible edge to its right may be extended 
to scanline j. Each visible edge that is not matched by a video partition is not allowed to 
extend to scanline j but rather is terminated on a previous scanline. At the conclusion of 
Ry scanlines, the visible edge encoder has encoded a set of visible edges that define the 
entire video image. 

The visible edge encoder maintains two major X-ordered tables. The first of these is 
the Current Visible Edge Table which holds visible edge data valid at the last processed 
scanline j' and which may possibly extend to the present scanline j. The Xordering of 
the Cunent Visible Edge Table is that of scanline j'. This X-ordering is maintained in 
a companion table that also maintains a listing of free locations in the Current Visible 
Edge Table. The size of these tables is 256 N- 1 where N is the number of memory banks 
in the CHARGE terminal. Each entry in the Current Visible Edge Table is of the form 
of Figure 65. Sg,Sj^, S last TYPE> and M are as defined in the video partition originating 



S|ast 
TYPE 

M 

X 

be 
PNTR 

Figure 65 - Current Visible 
Edge Table Entry Format 
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the visible edge* X and are the horizontal position and the brightness of the edge at 
scanline j*. PNTR is a pointer to tue location of the edge in the second table: the Current 
Edge Block Table. 

The Current Edge Block Table holds visible edge data for all valid edges in an edge 
block. Some of these edges may still be current and have entries in the Current Visible 
Edge Table, while other edges may have ceased to be valid and are flagged as "done." 
Each entry in the Current Edge Block Table is of the form of Figure 66. 



COLOR 



Figure 66- Current Edge 
Block Table Entry Format 



is the scanline on which the visible edge was first created, COLOR is the color to the 
right of the edge; is the horizontal position of the edge at Y = Y^; and b^ is the brightness 
of the edge at Y = Y^. The remaining parameters are filled in when the edge ceases to be 
"current** and is flagged as "done." H is the height of the edge such that Y = Y^ + H is the 
last scanline for which the visible edge is valid. M* and g* are representations of the edge's 
slope and brightness gradient and are given by 

M* - + [3.11] 

g* = [{bc-bt)2l^-P/H + .5] [3.12] 

where p is the number of high order O's in an 11-bit H. If the surface to the right of the 
visible edge is to be smooth-shaded (i.e., TYPE = 1 or 3), the value of Y^ is set negative 
flagging this condition. Companion to the Current Edge Block Table is a second table of 
equal lengcn (256 N 1) which maintains the X-ordering of the visible edge block, a listing 
of free locations in the Current Edge Block Table, and the "done'' list that flags those 
visible edges that have been completed and had values for H, M*, and g* calculated. The 
X-ordering of the Current Edge Block Table is from left to right (high X to low X) and 
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is such that for any scanline, those edges valid on that scanline are X-ordered between 
themselves. Tne X-order between groups of edges disjoint in Y-ranges is arbitrary. 

After a block of visible edges is complete, the visible edge encoder sends it to an 
output table in sequential left-to-right order preceded by a marker edge indicating the 
Y-range over which the block of edges is valid (see Figures 46-48). The size of the 
output table is the size of the CHARGE terminal's memory. 

The last table that the visible edge encoder uses is a sequential right-to-left X-ordered 
Scanline Output Table remembering the X-order of current visible edges as they are 
encountered on Scanline j. Each item in this table has three entries: 

(1) The location of the edge in the Current Visible Edge Table 

(2) A bit to say if the edge is "new" (i.e., if the edge has been created 
on Scanline j) 

(3) A pointer to the edge in the Current Edge Block Table that is the first 
edge to its right on scanline j. This table is used for updating the 
X-ordered pointers for the Cunent Visible Edge Table and the Current 
Edge Block Table after the matching processes have been completed. 

Figure 67 provides a functional flowchart of the operations on scanline Y = j\ 
If j = 1, the display is initialized clearing out all of the vi^jible edge tables and setting 
^last = 0 (where Yj^gt is the last scanline processed by the visible edge encoder) and 
^block ~ 1 (where Yy^^j^ is the initial valid scanline for which the current edge block 
is valid). Next, A Y = j - Yi^|. allowing positional updates of current visible edges 
(Xi = Xj + A Y-Mj). A test as to whether or not the current edge block should be 
ended and another begun is then made. This test is 

N-(256- [(Yb,ock-J)/40) 1)>Ncb + Nv [3.13] 
where the left side of the inequality represents the maximum number of non-marker 
edges in a visible edge block (see [3.3] ), and the right side represents the maximum 
number of potential visible edges in the current edge block at the end of proceosing 
for scanline Y = j. Here, Ncy is the number of visible edges presently in the Current 
Edge Block Table (both current edges and "done" edges) and Ny is the number of 
video partitions. Obviously, none of the visible edges may match any of the video parti- 
tions yielding Ncg + Ny as the maximum possible number of visible edges at the end 
of processing for scanline Y ^ j. Initial scanline processing is completed with a check 
for end of display (i.e., Y > Ry). If the end of display has been reached, the current 
edge block is completed and sent to the output table. Then the rest of the output 
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table is filled with null edge blocks that serve both to (a) fill up the remainder of the 
CHARGE terminal's memory and (b) provide null data for those Y values (greater 
than Ry) that represent vertical retrace time between interlaces. 

Once initial scanline processing has been completed, the major processing loop 
of matching the video partitions and the visible edges is entered. The procedure is to 
attempt to match each video partition (working in X order from right to left) with 
an unmatched visrble edge (scanning from right to left). The parameters used in matching 
are Sg, S^, S\^i, and TYPE. If no match is possible, a new visible edge is created. If the 
visible edge that is matched is not the first of the previously unmatched visible edges 
(X-ordered from right to left as they appeared on scanline j', the last scanline worked 
on by the visible edge encoder), all of the visible edges up to the matched visible edge 
are terminated on scanline j- 1 since any subsequent attempt to extend one or more 
of them would lead to an illegal ''crossover" condition (i.e., there would not exist a 
valid X-ordering of the edges in the Current Edge Block TabJe), The brightness and 
the horizontal position of the visible edge are updated to scanline j by setting 
be = bp and X = X + A Y-M. Of course, if the visible edge has just been created on 
Scanline j, then X = [Xp] • 

The CHARGE terminal uses M*, a less precise value of M, as it decodes the 
visible edges into an image. In Appendix C it is shown that for Rx < 2048 and 
Ry < 1024, a 13-bit M* (12 bits plus sign bit) can lead to a maximum deviation oi 
e - 0.5 units in X, Note, too, in Figure 68 that there is an error or up to 1 unit in 
X in the position of X^ between the original projected edge and the visible edge* The 
error is due ,to the necessity of starting the edge at the CHARGE terminal (and, hence, 
the visible edge encoder) at the discrete X = [Xp] rather than at X = Xp. The design 
of the CHARGE terminal is such that decoded edges may never cross each other; they 
may come as close as 1 unit in X, but they should maintain an average separation of 
two units in X between every three edges. 

Accordingly, after an edge has been matched ;nth a video partition and tentatively 
extended to scanline j, it must be compared in closeness to the last visible edge to its 
right to make certain that there is no danger of a crossover at the CHARGE terminal* 
Hence, after an edge has been matched, it enters the processing of Figure 67. Two 
registers are used to remember X information about previous edges on the scanline. 
XLST contains the X position of the last vL.ible edge to the right; XLTCH contains 
the X position of the last visible edge to the right that satisfied the requirement of an 
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Figure 68 - Projected Edge; Edge at Visible Edge Encoder; and 
Edge at CHARGE Terminal 



average separation of two units in X for the last three edges. If (XLST-XLTCH) < 3 
where XLST-XLTCH, ILTCH points to the edge at XLST which has been tentatively 
extended to scanline j but which may be deleted if the present edge, positioned at 
X = Xc, is "toe close." 

Figure 69 demonstrates the normal mode of operation on scanline Y = j. The la^t 
line worked on was Y = Yjagt. Seven scanline partitions on line Yj^^t generated seven 
visible edges in the Current Visible Edge Table. These edges have the order: a, b, c, d, e, 
f, g. If no other visible edges are present in the Current Block Table, the entries in Current 
Block Table have the order: g, f, e, d, c, b, a. Note th.'^t a visible edge/scanline partition 
always is present at X = 1 and X = Rx, due to the 1:0 and 0:1 projected edges present in 
every image. 

As the processing of scanline Y = j begins, the first partition immediately matches 
current visible edge "a." The next two video scanline partitions immediately match current 
visible edges "b'' and "c/' As each visible edge is matched, it is removed from the X-ordered 
list of current visible edges yet to be matched. The fourth video scanline partition is then 
compared with the next visible edge, "d," and it is found that they do not match, however. 
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Figure 59 - Visible Edge Encoding on Scanline Y=j 

it is found that partiti n "4" matches visible edge "e." Since edge "e" is not the first 
visible edge in the X-ordered current visible edge list, it is deemed that those edges with 
lesser X-values than edge "e" (here, only edge "d") are no longer valid and should be 
terminated on scanline Yjj^g^. 

Next, partition "5" is encoded. If partition "5" is tagged as "new" (i.e., it is the 
top of a projected edge), a new visible edge is immediately created. If, however, it is not 
tagged as "new,'' then partition "5" is compared with all of the remaining current visible 
edges in search of a match before a new visible edge is created. 

Finally, partitions "6" and "7** match current visible? edges "f " and "g." The 
Current Visible Edge Table now has the X-order: a, b, c, e, h, f, g. The Current Edge 
Block Table has the X-order: g, f, h, e, d, c, b, a. Figure 70 demonstrates the X-ordering 
of a block containing edges for a cube. 
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VERTICAL SKIP.AHEAD 



There is, in general, a great deal of coherence between successive scanlines on which 
no new projected edges become valid or no projected edges have ceased to be valid. This 
is due to the fact tliat a surface may as well be thought of as generating vertical Y-ordered 
partitions over which it is valid for a given horizontal X value, rather than horizontal 
X-ordered partitions over which it is valid for a given vertical Y value. Simply stated. 
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Figure 70 - Order of Edges for Block Depicting 
a Cube. 



visible surfaces tend to remain visible unless the edges defining them end or unless new 
edges defining new visible surfaces begin. 

Each scanline processing step requires the following operations: 

(1) Adding new projected edges to the set of current projected edges. 

(2) Updating the current projected edge set for the next scanline. This 
processing includes deletion of edges that cease to be valid (i.e., H goes 
negative upon update). 

(3) Generating actual video scanline partitions. 

(4) Possible readout of a completed visible edge block. 

(5) Updating current visible edge set. 

Furthermore, it is known that at the completion of a scanline processing step, the 
actual X-ordered video partitions are matched exactly by the current set of visible edges 
if no deletions of visible edges occurred due to the "closeness'' comparison test. Hence, 
if the order of the actual video partitioning of the scanline remains constant over a skip- 
ahead attempt of A Y scanlines over which the set of valid projected edges is constant, 
then the skip-ahead attempt is known to be successful. 
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If no edges start or end on a given scanline, the above set of operations reduces to: 

(1) Updating the corrent projected edge set by updating Dq and X and 
maintenance of the X-ordered edge list. 

(2) Generating the actual video scanline partition. 

(3) Updating the current visible edge set. 

Vertical skip-ahead attempts require the following processing: 

(1) Saving the current projected edge set and its X-ordered list. 

(2) Saving the actual video scanline position. 

(3) Updating the current visible edge set by updating Dq and X and generating 
an X-ordered edge list. 

(4) Generating the actual video scanline partitions. 

(5) Comparing this new scanline partition with the saved scanline partition. 

(6) If the skip-ahead attempt is successful, the parameters of the projected 
edge set and the X-ordering of the projected edge set must be replaced by 
those that are valid at the skip-ahead Y value. 

If multiple skip-ahead attempts are made v/ithin a range of scanlines, operations (3), 
(4), and (5) are repeated with each attempt. 

Generally speaking, the average amount of time required to process A Y scanlines 
without any vertical skip-ahead is given by 

To = AY-E[to(n)] [3.14] 

where 

n : Number of edges in the current projected edge set. 
to(n) : Random process mapping n into an estimate of processing time required 
to process n projected edges into a^^tual video scanline partitions and then 
into X-ordered visible edges. 
E('] : Expectation operator. 

If a vertical skip-ahead is successful over the same A Y scanlines, the average amount 
of time required is given by 

Ti=Ts(n) + E[ti(n)] [3.15] 

where Ts(n) is a setup time and ti(n) is a random processing time (E[ti(n)] > E [to(n)] ). 
Let two hypotheses be: 

Hq : Vertical skip-ahead will be successful with probability P. 

Hi : Vertical skip-ahead will not be successful with probability (l-P). 
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The cost matrix for the vertical skip-ahead decision problem is of the form 

.^10 ciij 

where cy is the cost of choosing Hj when Hj is actually true* Here, 

"to to 



[3-161 



[3-171 



The average risk is given by 



C= PlcQoPr (choose Ho, Ho true) +cioPr (choose H^, Ho true)) [3-18] 
+ (1-P) IcQiPr (choose Ho, H^ true) + ciiPr (choose Hi,Hi true)] 
The problem of minimizing [3,18] is quite complex- Two approaches have been 
considered. 

APPROACH!: If the joint probability density functions (p.d.f-) Po (n, Ay) and 
Pi (n, Ay) for the two hypotheses Hg and Hj^ were known a pnori, it would be 
possible to form the likelihood ratio 

A(„,Av).^li=44 (3.19] 
Po (n. Ay) 

which would be compared to the threshold value 



^ (1"P) (coi"Cii) (1-P) To"Ti 



[3-201 



The strategy would then use the Bayes decision rule 

A< Ao , chose Ho 

A>Ao , chose Hi [3,21] 
A = Aq , chose either Ho or H^ 

Hence, to implement APPROACH I it is necessary to analyze a vast number of 
images, gathering data on values of (n. Ay) for which vertical skip-ahead was and was not 
successful. From these data, an empirical decision surface could be defined, (see Figure 71). 

APPROACH I has not been implemented because it was feared that any particular 
sample of test images would not be sufficiently varied as to converge upon the true Bayes 
risk; that is, any particular set of images would unintentionally weight the decision surface 
in one way or another. 
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' Choose Hq 




Choose Hq or H-j 



Figure 71 - Empirical Skip-Ahead Decision Rule. 



APPROACH II: The second skip-ahead approach seeks to make use of certain known 
properties of projected atoms and c acts via a sophisticated, but suboptimum, skip-ahead 
strategy. These known properties are: 

(1) If n = 2, only the two border edges at X = 1 and X = Rx arn present. 
Hence, if n = Z, skip-ahead will always succeed. 

(2) For n large, Tj » Tq/A Y; that is, the time required in a skip-ahead 
attempt becomes much larger than the average time per scanline for normal processing. 
Also, if n becomes large, the probability of placing n projected edges on a scanline 
without changing the order of visible edges becomes small. Hence, if n > fl^* a threshold, 
vertical skip-ahead attempts are prohibited. 

(3) Since Tj > Tq, for A Y < Oy, a threshold, it does not really pay to use 
vertical skip-ahead. 

(4) Actual video scanline partitions change from scanline to scanline because 
either (a) the order of visible surfaces changes or (b) projected edges come **too close." The 
former case is roughly related to the number of distinct atoms in the region A y (and 
hence is related linearly to n); such crossover cases occur with equal probability 
throughout A y. The latter case occurs at the vertices of atom facets (i.e., scanlines on 
which projected atom edges are started or are terminated). Thus, a reasonable probability 
density function for a successful vertical skip-ahead might look like Figure 72A, while 
possible p.d.f .'s for vertical skip-aheads that fail are like Figure 72B-E. That is, if a 
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Figure 72. Probability Density Functions for Vertical Skip-Ahead. 
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s ap-ahead attempt of A Y = Yj scanlines should fail, the most likely place for it to fail 
(i.e., pi(n, Ay)/po(n, Ay) = max) is at one or both end points in Y^. 
Hence, APPROACH II uses the following algorithm: 

(i) If n = 2, always perform a skip-ahead. 

(ii) If n > 0^ never attempt a skip-ahead. 

(iii) If A Y < d y, do not attempt a skip-ahead. 

(iv) Otherwise, attempt a skip-ahead of the full A y = Yj scanlines. 

(v) If an attempted skip-ahead should fail, proceed according to the 
functional flow of Figure Ti. This algorithm, having first failed to accomplish the maximum 
skip-ahead (Figure 72A), successfivly attempts to decide whether pi (n, A y) is of the type 
shown in Figure 72B, 72C, 72D, or 72F. That is, having failed at A Y = Yj, it next attempts 
at A Y = Yi/2 should fail, the p.d.f . is probably either type C, D, or E and the next attempt 
is at A Y = 0 y. The process continues each step remembering what has been tried before 
until the piocess is complete. At each stage of the process, rule (iii) applies. 
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Appendix A 
EXAMPLE OF A NON-PLANAR SURFACE 

Consider the right cylinder of radius R and height H in Figure A-1. 
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Figure A-l - Right Cylinder' 



The points on the curved Surface S, satisfy 



0 0 0 H H 
f{x) = 0 = + X3 - R^ "^^^2^2* 



The gradient vector at any of these points on S is given by 

NM = Af"^(x) = (2x^,0,2x3)"'" 



[A-1] 



[A.2] 



The problem is to represent the closed surface S with a set of n>2 planar surfaces. The 
approximation error at any point x on the original surface S will be dei, led as the dis- 
tance between X and x, the closest point to x on the approximating surfaces jSjj . 
That is 



e{x) = minjll X- Xjllj x. e Sj | . . n 



[A.3] 



Since the gradient vector N(x) for S happens to be independent of X2, any surface 
Sj representing a portion of S should also have a gradient vector NjCSj) which has a zero 

A A 

X2 component to minimize e(x) over Sj, that is, any surface Sj may be extended in the 
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X2 direction without incremental error (i.e., '^thout adding to the norm II x- xjll). This 
means that the most efficient set of n planar surfaces is a set of n rectangles as in 
Figure A-2. The rotational separation between adjacent facets is fl=2ff/n=ff/4. Also note 
that rather than just the end points, each of the vertical edges in Figure A-2 lies entirely 
in S. This is an additional benefit of having N(x) independent of X2 (i.e., curvatic:i in 
only 2 of the 3 dimensions). 
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Figure - Right Cylinder Represented by 
n = 8 Planar Surfaces 



For the right cylinders of Figures A-l and A.2, the error for any one facet Sj (at 
X2 = constant) is as shown in Figure A-3. x is a point on S; x is the closest point on S; 
and x^ and X2 are points on the vertical edges of S. S has been arbitrarily located such 
that it is bisected by the X3 = 0 plane. The error is given by 

e(0) = llx"5ll 

= llRcos0- Rcos0/2ll -e/2<<l> <0/2 
el0) = R Icos0 - cosO/2 I - 0/2 < 0 < 0/2 
In terms of n rectangular facets 

e(0) R{cos0 - cos7r/n) -7r/n <0 <7r/n [A-4] 
and the maximum error over any one facet is given by 

W ^ ^ '^(l - cos7r/n) [A.5] 
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Figure A-3 - e(0) = e{x) = II X - 
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The normalized error of the approximation of S by |Sj| is given by ejyj = e/R. 

= cos 0 - cos ff/n - tr/n < 0 < tr/n [A-6] 
Gm = 1 - cos tr/n [A-7] 



'max 



Figure A4 demonstrates the reduction of ew as the number of approximating 

max 

surfaces is increased. It should be noted that after an initial rapid improvement in the 
approximation, the marginal improvement in the approximation becomes quite small. 

The brightness at any point x on Surface S due to a sunlight source at infinity 
(i.e., parallel light rays of constant intensity through R^) is given by 



° \llN(x)il llxJI 



[A-8] 



where Bg is a scalar measure of the sunlight intensity, Xg is the direction of the sunlight 
through r3, N(x) is the normal vector to surface S at point x, and f(a) is the brightness 
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Figure A'4 - Reduction of ^nmax by Increasing the 
Number of Cylinder Facets 
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reflectivity function for surface S. A typical brightneso reflectivity function for a 
nonshiny surface is given by 

+ A2 + cos a 



f(a) =< 



A . A . . lal<7r/2 

A A A [A-9] 

Al + A2 + A2 cosa TT 
if !al>rr 



A^ + A2 + 1 2 



where and A2 represent (as fractions of the direct light reflectivity) the effects of 
diffuse light in R^. Aj^ is an overall ambient light that affects the entire object; that is, 
it adds light to all surfaces equally whether or not they face toward the sun. A2 repre- 
se**ls the effects of reflected light from surrounding objects and provides the smooth- 
shading of the dark -side object surfaces (i.e., those surfaces that face away from the sun). 
Hence, the reflectivity function, f(a), falls off as the cosine of a from a peak value of 
1 at a = 0 to a value proportional to A^ + A2 at a = 7r/2. As a increases from 7r/2 to tt, 
f(a) falls off as the cosine of a to a baseline value proportional to A^. 

If Ai = 0, those surfaces pointing directly away from the sun are black. If A2 = 0, 
all surfaces facing away from the sun have equal brightness. If A^ = A2 = 0, the effect 
to the eye is as if the object were in outer space with bright sunllL buxfaces, and black 
surfaces which face away from the sun. Typical values for Aj and A2 on a clear, bright 
day are Aj = .05 and A2 = .3. Figure A-5 demonstrates the behavior of [A-9] . 

Again consider the cur/ed surface S of the right cylinder of Figure A-1 and the 
set of approximating surfaces jSjj of the approximate right cylinder of Figure A-2. 
Due to symmetry about the X2 axis, the brightness shading of both S and jSj | is 
independent of the X2 component of x^. Hence, the simpler, 2-dimensional case of 
Figu^'e A*6 may be considered. Xg and N'(x) are the (normalized) projections of N(x) 
and Xg onto a constant X2 plane. The angle 0 between Xg andl^'(x) is given by 
cos" 1 (N'T(x) • x'g). The brightness at any point x on S is given by [A-8] and is 
merely a scaled plot of f(0). However, the precise brightness for jSij is known only at 
the n points (equally spaced, but randomly phased W.R.T. (j>) in 0t [0,27r ] . The brightness 
between pairs of these n points is given by either the linear rule 

B(X) = B(i) +X (B(i+1) - B(i)) \ e [OJ] [A-IO] 

or the logarithmic rule 



B(\) = B(i) 



/Bji+r 
VB(i) 



J xeioj] [A-llJ 
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The actual brightness, B(0), and two approximating linear brightness functions, 
8^(0 ) and 82(0), phased the maximum distance of 0 - n/n - n/S radians apart, are 
plotted in Figure A-7. 

It should be noted that both Bj[(0) and 82(0) use the same brightness approximation 
rule in shading |Si| over n'^S points in the angular range [0,27t]. It is simply the chance 
orientation of the 8 facets W.R.T. the sun vector that causes Bi(0) to be a closer approxi- 
mation than 82(0). Not only does §2(0) fail to reflect the true brightness extremes of 
8(0), but it also erroneously indicates regions of constant brightness at 0 = ± krr, 
k = 0, 1, 2, . • . . Hence, 82(0) is the worst case linear estimate for B(0) and therefore 
should be used in any brightness error evaluations. 

to 
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The maximum errors between B{(t>) and 62(0) occur at ^ = 0 and ^ = ± 7r/2. The 
brightness error at ^ = 0 is given by 

65(0,0) = 85- Bscos-" 
n 

eb(0,n) = Bs (1 - cos ^ ) [A-121 

which has the same form as [ A-7] , the error in the surface approximation. The brightness 
error at ^ = - n /2 is maximum when A2 = 0 resulting in constant brightness across the dark 
side of the object, Figure A-8. 




Figure AS-- Brightness Error eb(-7r/n, n,0i) 
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The bri^itness error at 0 « - 7r/2 for a phase angle of 0i radians is given by 

^,..„,„,,„ = |il.(i.„) ,A.13, 

where 0i e [0, 7r/2n] . Since eb is zero at 0i = 0 and 0i = Tr/2n, the maximum value of 
eb (-^/2, n, 0i) can be found from 

— = 0 = -sin0i + (7r/n-0i)cos0i [A-14] 

or 

tan 01 +01 = TT/n 

For n sufficiently large tan 0i = 0i ^ 7r/2n. For example, if n = 8, 01^ tr /2(8.1). Thus 

, . V . s«n Tr/2n ^ ^ r A i Ki 

eu (- 7r/n, n) = n > 8 (A-ISJ 

2 

Table A-1 lists values for eb(o,n) and eb (- 7r/2, n) for various values of n. As expected 
eb(- TT/n, n) > eb(o^) since the Taylor series expansion for eb(- n/n, n) is dominated by 
(7r/n)2/2!, while the Taylor series expansion for eb(o,n) is dominated by the smaller 
(^/n)3/3!. 



n 


eb(-W/n,n)/Bg 


eb{0.n)/Bs 


8 


.19134 


.07612 


9 


.17101 


.06031 


10 


.15451 


.04894 


20 


.07821 


.01231 



Table A-1. Brightness Approximation 
Maximum Errors 



The human eye is also quite sensitive to the rate of brightness change. The derivative 
of the brightness functions of Figure A-7 is shown in Figure A-9. The derivative of the 
true bri^tness function B'(0) has discontinuities only at 0 = ± 7r/2 in 0€ (- n, n] ; this 
discontinuity occurs on the locus of points separating the sunlit side of S from the side 
of S illuminated only by diffuse light. Hence, that discontinuity serves a very useful 
purpose to the eye* 

It can be noted from Figure A-9 that B'j (0) appoximates B'(0) more closely than 
B'2(0V Both B'i(0) and B'2(0) converge to B'(0) as n becomes large. 
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Figure A^Q-- Brightness Function Derivatives 



The same general conclusions occur if the logarithmic rather than the linear smooth* 
shading rule is utilized. 

LinearRule: B(X ) = Bq +X (Bf- Bq) Xe(0,l) [A-16] 



er|c 



Logarithmic Rule: B(X) = BqIBj/Bo)^ X€[0,ll 
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[A-17] 
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The deviation of the two approximations is given by 

D(X) = B(X)-B(X) 

= Bo +X (Bf - Bq) - (Bf/Bo) ^ [A-IS] 

and the normalized deviation is given by 

Dn{X) = — = 1 + X(^-^ - l) - Bf/B/ {A.19] 

*rhe value of X for which X ) is maximum can be obtained by setting 9Dn(X )/9 X =0 
and is given by 



£n (Bf/Bp-I) -£n (gn {B^/Bq)) 

^max ~ 



Bn (Bf/Bo) 



[A-20] 



and so X : V4 1 as (Bf/Bo) : l"*^ as shown in Figure A-10, Figure A-11 demonstrates 
the deviation of the approximations for (Bf/Bo) = 2 and (Bf/Bp) = 10. 
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Appendix B 
LOGARITHMIC BRIGHTNESS PARAMETERS 



The problem here is to logarithmically encode a brightness h€ [Bj^in^^maxl ^^^^ ^ 
m-bit binary number. Here, B^ax > ^xnin > 0- % ^^e logarithmic rule 

B - B^j, (B^3,/B^;,y [B.l] 

for some XG[0,1]. Then 

B 



Bmin \Bmin / 



^ xeio.il [B.2] 



and taking logarithms on both sides of [B.2] yields 



^min / 



From a computational point of view 

12"" - 1) 

b= g^TB-mlxX fp— \ ' bei0,2'^ - 11 



'2"" " / n \ 



(2"" - 1) 
[l!nB^ax-fi"Bminl 
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[B.3] 



and so 

(2'^-1)logA. b bGlO,2'^-11 [B.4] 



b=;r-7r — r f^n B- finB^J bel0,2'^-ll [B.5] 



or simply 

b-KifinB-K2 bel0,2""-ll [B.61 

If Bmin = 1. [B.5] reduces to 

12""- 1) 

b==— finB belO,2'^-H [B.7] 

Typical values are B^in = 1, B^ax = 149. m = 7 and so 

127 

b=— il-enB be 10,1271 [B.8] 
tn 149 



Appendix C 

ERRORS IN VISIBLE EDGE POSITION AND BRIGHTNESS AT TERM 



The CHARGE image generator uses 4-byte floating point representation for M-the 
dope of the visible edge, and g-the gradient of brightness along the edge. The CHARGE 
terminal is sent lower resolution representations of the slope and brightness gradient 

^""'^ M-.U [m-2"-.,5] 2-" tC-U 

[,.2 "-P.. 5] 2"-" 
where the operator tr [•] indicates truncation to the next lower or equal integer. In both 
[C-1] and IC-2], p the number of high-order zeroes in an 11-bit fixed-point 
representation for the edge height, H. Since in general M M* and g g*, there is 
generally a discrepancy between the bri^tness and horizontal position at the CHARGE 
image generator and at the CHARGE terminal. In particular, since edges must not cross 
each other at the CHARGE terminal, allowances for the slope error must be made at 
the image generator. If two edges approach each other to within the possible error 
amount, one must be terminated. 

The horizontal error of any one edge at a distance Ay from the top of the edge 
is given by 

ej^(Ay) = + M-Ay - (x^ + M*-Ay) 

e^^lAy) = (M- M*)Ay 1^-3] 
So lex(Ay)l attains it*s maximum value when Ay = H. 
Hence 



e^< Ax 



and so 



le,l<l2P-^2^l ICS] 

1 1"* D 

If p = # high-order zeroes in an 11-bit H, then H < 2 ^, and so 

Ie,l<0.5 lC-6] 
where 1 = 1 horizontal discrete x- resolution interval (1/1600 of TV width) 



Hence, since ex^( - 0.5, 0.5), two edges should not be allowed to come within one 
x-unit of each other at the CHARGE image generator to avoid the danger of crossover 
at the CHARGE terminal. If a minimum separatioi. of 1 horizontal resolution unit is 
desired at the CHARGE terminal, the closeness criterion at the CHARGE image generator 
must be 2 units. 

A similar analysis for the brightuess gradient yields a maximum brightness error of 

leb 1< 0.5 [C-7] 
where 1 = 1 brightness resolution interval (1/128 of brightness range). 




Appendix D 
CHARGE COMMANDS 



Each CHARGE command falls into one of the following categories: 

1. Atom commands— used by an author to create different 3-dimensional 
shapes (atoms). 

2. Object commands— used by an author to define complex objects in 
terms of atoms and previously defined objects. 

3. Viewing commands— used by an author to define parameters to be used 
by the image generator when displaying a collection of objects. 

4. Image commands —used by an author to store and retrieve the results of 
an image generation process. 

Commands may be entered in a free-form format with blanks separating the 
various components of each command and the character *7" serving as a delimiter 
between successive commands. 



ATOM COMMANDS: 



(1) 
(2) 
(3) 
(4) 
(5) 
(6) 
(7) 



AT/ 

C Atom/ 



Invoke atom command processor. 



Create Atom. 



M Atom/ 
XAtom/ 
LA/ 
U/ 



C Atom 1 



Atom 2/ 



Create Atom 1 using Atom 2. 
Modify Atom. 

Scratch Atom from uncompiled Atom file. 
List atoms in uncompiled Atom file. 



Compile atom currently being modified, and 
put it into compiled Atom file. 



(8) 



A (cm) . . 



/ 



Add atom components: Point:P(#,x,y,z) 
Surface: S(#,Pi,P2,P3) 
Brightness Normal: B(# ,81,82 . • .,810) 
Edge: E(#, Type, Pi,P2,Sl,Sr,Bi,B2,B3,B4) 



(9) 
(10) 

(11) 



D (cm) . . 

LP/ 

LB/ 



/ 



Delete atom components. 

List all atom points. 

List all atom brightness normals. 




ATOM COMMANDS (Cont.): 



(12) 


LSI 




TiiQf" 2)11 ntntn ciirfippc 

JLJlOu CUl CXLV^lil oUXlclV^V'O* 


(13) 


LEI 




T.icf nil ftfniYi orlcToc 


(14) 


LP ni - 


no/ 


T.icf afnm nnitifc n-i fn no 


(15) 


LB ni - 


no/ 


Til^f nfnm hHcllfnPQQ nnrmolc n-i no 
JLJlOv c&u^lll Lll Jgli LliCoo llUxlllulo 111 \AJ I1V« 




JuO 111 




List atom surfaces n^ to n2. 


(17) 


LEni- 




List atom edges n^ to n2. 


(18) 


AL/ 




Invoke object command processor. 


(19) 


IL/ 




Invoke image library command processor 



0BJEC7' COMMANDS: 

(1) AL/ 

(2) C Object/ 

(3) M [Object]*/ 

(4) S [Object]/ 

(5) X [Object]/ 

(6) A Object [tr]/ 



(7) D[ object [.#]],,,/ 

(8) I Object [.#] #[-#] [tr] 

(9) R Object [.#]#[-#] [tr], 

(10) O [Object]/ 

(11) T [Object 1 [Object 2] 1/ 

(12) N Object 1 Object 2/ 



Invoke object command processor. 
Create object. 
Modify object. 
Show object. 
Scratch object. 

Add sub-object with transformations: 
m(x,y,z): Move in x,y,z 
S(x,y,z): Scale in x,y,z> 0 
X(x,y,z,Ox): Rotate in X about (x,y,z) 
Y(x,y,z,^;y): Rotate in Y about (x,y,z) 
Z(x,y,z,02): Rotate in Z about (x,y,z) 
OR(P): Optical reflectivity Pe [0,1] 
C(b,g,r): Color if not yet colored. 
RC(b,g,r): Color 
MIX: Mirror image in X. 
MIY: Mirror image in Y. 
MIZ: Mirror image in Z. 

Delete sub-object. 
./ Insert transformations. 
./ Replace transformations. 

List object structure. 

List transformations. 

Name Object 1 Object 2. 



♦An argument in brackets denotes an optimal argument. 
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OBJECTXmiMAimSJSlont\ 

(13) L [#] / List object names. 

(14) LA/ List atom names. 

(15) UC Object/ Show use count of object. 

(16) AT/ Invoke atom command processor. 

(17) IL/ Invoke image library command processor. 

(18) ZZ/ Drop CHARGE system. 



VIEWING COMMANDS: (processed by object command processor) 

(1) V [par] . . ./ Set viewing parameters. 

(2) VR [par] . . ./ Reset viewing parameters to their default 

values and then set parameters: 
W(x,y,z): Viewing window locp.tion. 
E(x,y,z): Eye location with respect to win Jew. 
0(x,y,z): Point observed. 
ws(x,y): Window size in ' .^ters. 
MAG(M): Magnificatior M>0. 
FS(F): F^stop for observer F€ [0,1] . 
BR(B): Maximum brightness at terminal. 
RES(Rx,Ry): Resolution of tenriinal. 
SUN(x,y,2): Sunlight direction. 
SB(.5): Fraction of B for sun. 
BC (b,g,r): Background default color. 
BB(b): Background default brightness. 
AB(ai,a2): Ambient brightness factors. 
T(t): Terminal type: 
1 : normal 

2: no smooth-shading 



IMAGE LIBRAE Y COMMANDS: 

(1) IL/ 

(2) SU Image/ 

(3) SH Image/ 

(4) LI/ 



Invoke image library command processor. 
Save Image. 
Show Image. 
List image names. 
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IMAGE LIBRAE Y COMMANDS (Cont.) 



(5) 


SN [#1 . . ./ 


Show imaEfefs^ hv nnmhpr 


(6) 


DI Image . . ./ 




(7) 


RI Image 1 Image 2/ 


RGnarriG Imai?G 1 Ima{?G 2 


(8) 


CP Packet [Image] . . ,/ 


CrGat*^ imaiTfe nackGt. 


(9) 


DP Packet/ 


Dplpfp narkpt 


(10) 


LP Packet/ 


List oacket contents 


(11) 


SP Packet/ 


Send oacket. 


(12) 


RP Packet 1 Parkpt 2/ 




(13) 


RS/ 


Resend last image. 


(14) 


AT/ 


Invoke Atom command processor. 


(15) 


AL/ 


Invoke Object command processor 
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